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ABSTRACT 



Initial yield surfaces of unidirectional fiber reinforced 
composites subjected to combined longitudinal, transverse 
normal, and longitudinal shear loads are calculated. The 
composite is composed of filaments arranged in a periodic 
square array embedded in a matrix material. The constituent 
materials are assumed to be homogeneous and isotropic, while 
the composite is assumed to be macroscopically homogeneous 
and transversely isotropic. The finite element method, 
using linear strain triangles, is employed to calculate the 
stresses throughout the composite. Von Mises yield criterion 
is used to calculate the elastic limit of the local micro- 
scopic combined stresses in the composite. A parametric 
evaluation is carried out by changing individual constituent 
properties and evaluating the effect of this variation on 
the composite properties. Yield surfaces for several 
functional composites are included, as well as a computer 
program, in FORTRAN IV language, for the calculation of a 
yield surface. 
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I. 



INTRODUCTION 



A. COJ-IPOSITE I4ATERIALS 

The requirements of modern technology, especially 
that of the aerospace industry, have necessitated the 
development of materials possessing much higher strength 
to weight ratios than common structural materials. This 
result has been achieved through the introduction of com- 
posite materials, that is two or more constituent materials 
combined in such a way that the gross properties of the 
composite are favorable combinations of those of the 
constituents. This indicates that the desirable properties 
such as strength/weight and stiff ness/v;eight ratios are 
increased while less desirable properties such as .low 
ductility are minimized. 

A specific type of a two constituent composite is a 
unidirectional fiber reinforced composite. This consists 
of continuous lengths of a stiff reinforcing fiber imbedded 
in a ductile matrix material. All the reinforcing fibers 
are oriented in the same direction within a single layer of 
material. Generally several layers, or laminae, are bonded 
one on top of another to form a high strength laminated 
composite. 

Each of the lamina may be considered as an anisotropic, 
homogeneous material. The transverse properties, i.e., 
those in a direction' perpendicular to the filament orienta- 
tion, are weaker than those in the longitudinal (fiber) 
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direction for a single lamina. To counteract this weakness 
in a laminated composite, adjacent lamina are arranged such 
that they have different filament orientation. The funda- 
mental unit is therefore the single layer of unidirectionally 
reinforced materials. It is this type of basic layer that 
is analyzed in this investigation. 

The physical dimensions of such a lamina can vary 
considerably, although the largest is generally but a fraction 
of an inch thick. The fibers vary in diameter from .0001" 
for glass and .0003" for graphite to .004" for boron. Lamina 
thickness depends on the materials used and the intended 
purpose. An individual lamina can be produced in nearly any 
desired thickness. However standard laminae are becoming 
popular. For example, most graphite comes in a standard .01" 
thickness. Allowing for the presence of the matrix material 
such a lamina may contain in excess of 25 fibers across its 
thickness . 

B. DEFINITIONS 

The terminology of composite materials is often confusing 
due to a lack of precision in defining specific terms. For 
example, at times it is unclear whether a term such as tensile 
strength refers to the ultimate strength of a material, i.e., 
its failure point, or to the point of initial yield. In 
composite materials it is necessary to be especially clear 
by noting whether a term refers to a composite property or 
to a property of a constituent material. In order to avoid 
such difficulties the definitions of this section will be 
adhered to in this study. 
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1 . Macroscopic Quantities 



These are quantities of the composite material. 

They are calculated by assuming that the composite consists 
of anisotropic, homogeneous layers. 

2. Microscopic Quantities 

These quantities are determined on a point by 
point basis throughout the composite. The different 
mechanical properties of the matrix and filament materials 
are considered in the calculation of these quantities. To 
see the relationship between macroscopic and microscopic 
quantities consider Fig. (1). Figure (la) represents a 
body with microstress, in Fig. (lb) the macrostress 

S in the x-direction obtained from the microstress field 

X 

of Fig. (la) is 

a 



S 

X 







(i-i) 



o 

S is the macrostress for the face x = b, and represents 
the average value of the microstresses. 

3 . Initial Yield Strength 

That combination of macroscopic stresses at which 
a point in the composite first reaches its elastic limit, 
i.e., no plastic strain has occurred in the composite. 

This study is restricted to plane macrostress 
fields with components S„, S^, and S„„, where the z axis 
is along the length of the fiber. Specific cases of 
initial yield strength associated with this investigation 
are 
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y 




(b/Y) 



a) Microstress Field a (b,y) 



y 




b) Macrostress S 

Figure (1) . Relationship of tiicrostresses and Macrostresses 



a) Longitudinal Yield Strength, when the only 
macroscopic stress is in the axial direction, i.e., 

7^ 0, S = S = 0. 

Z ' ’ X X2 

b) Transverse Yield Strength, when the only macro- 
scopic stress is in the transverse normal direction, i,e., 

S 7^ 0, S^ = S ^ = 0. 

X ' ' 2 XZ 

c) Longitudinal Shear Yield Strength, when the 
only macroscopic stress is in the longitudinal shear 
direction, i.e., S^^ 7 ^ 0, S^ = S^ = 0. 
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4 . Second Deviated Stress Invariant, 

The second deviated stress invariant at a point 
is defined as [l] 




(Oy - + (02 - 




( 1 - 2 ) 



This quantity can be related to the yield strength of a 
material. 

5. Effective Microstress ^ a 



6 . Specific Strength 

The ratio of a particular initial strength term, 
as defined above, to the specific density of the composite. 
The strength considered must be specified, for example, 
specific longitudinal strength. This may be applied to 
composite materials as well as single phase materials. 

7 . Specific Modulus 



to the specific density of the material. In the case of 
anisotropic materials it must be specified which modulus 
is being considered, such as the specific transverse modulus 



The effective microstress, a, is proportional to 
the square root of the second deviated stress invariant, [l] 




+ 6 




(1-3) 



The ratio of the Young's modulus of a material 
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8 . Ultimate Strength 



That point in macroscopic stress space at which 
either a matrix element or a filament element actually 
fails, i.e., breaks , cracks , or separates. Due to the 
brittle nature of filament materials the initial yield 
strength and ultimate strength usually coincide for filament 
materials . 

9 . Initial Yield Surface 

In a general three-dimensional system the stress 
state is defined by six stress components, three normal 
stresses, S^, S^, and S^,and three shearing stresses, 

and Sy.^. Thus an initial yield surface may be repre- 
sented in macrostress space. Any macroscopic stress state 
within the surface means that the effective microstress at 
every point in the composite is less than that required for 
yield. Any stress state lying on the surface means that 
some point in the material has an effective microstress 
equal to the elastic limit. For brittle materials without 
a plastic region loading states outside this surface are not 
possible. For ductile materials stress states outside the 
initial yield surface (within the limits of the ultimate 
strength of the material) will result in plastic deformation. 

The three-dimensional yield surface considered in 
this study is a specific case of the general initial yield 
surface described above. Macroscopic stresses, Sy, and 

Sy .2 set equal to zero in the calculation of the yield 

surface. The resultant surface represents the yield surface 
due to plane macrostress loading conditions. 
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10. Stress Concentration Factor 



A stress concentration factor, denoted as SCF, 
is generally defined as the maximum microstress in the 
composite to the applied macrostress. In this study three 
particular cases of stress concentration factors are 
investigated. 

a) SCFX. Transverse loading stress concentration 
factor, associated with macrostress state (S^^, S^, = 



( , 0 , 0 ) . 



SCFX = 



(^x) 



max 



X 



(1-4) 



b) SCFZ. Longitudinal loading stress concentration 
factor, associated with macrostress state (S^^, S^, = 



'xz' 



(0, S^, 0) . 



(Cfz) 



SCFZ = 



max 



(1-5) 



c) SCFXZ. Longitudinal shear loading stress 
concentration factor, associated with macrostress state 



(^x' ^z' ^xz) = (O' 0, S^,). 



(^xz) 



SCFXZ = 



max 



( 1 - 6 ) 



xz 



11 . Transverse Stiffness, E, p 

The transverse stiffness of a composite is the 
macroscopic modulus of elasticity in a direction perpen- 
dicular to the fiber direction. 
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12. Longitudinal Stiffness, E 

L 

The 'longitudinal stiffness of a composite is the 
macroscopic modulus of elasticity in the fiber direction. 

13. Longitudinal Poisson's Ratio, 

” _ - - Jj 

The longitudinal Poisson's ratio is the measure of 
the contraction in the transverse direction due to tension 
in the longitudinal direction, i.e., , for the case 



(Sx, Sy, S^) = (0, 0, S^) . 

14. Transverse Poisson's Ratio, v 

£ .p 

The transverse Poisson's ratio is the measure of 
the contraction in the transverse direction due to tension 
in the same plane, that is = e^/e^, for the case 



(Sx, Sy, S^)= (Sx, 0, 0) . 

15 . Shear Modulus , G 

The composite shear modulus is the modulus of 
elasticity in shear due to longitudinal shear loading. As 
the transverse composite shear modulus is a dependent 
quantity any reference to composite shear modulus in this 
study will mean the longitudinal shear modulus unless 
specifically stated to the contrary. 

16 . Filament Volume Fraction, v ^ 

The filament volume fraction is the ratio of 

filament material volume to total material volume, 

y = Filament Volume 
^ Total Volume (I”7) 



Note that 0 £ £ 1.0. 
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17. Matrix Volume Fraction, v 
— m 

The matrix volume fraction is the ratio of the 
matrix material volume to the total material volume, 

y _ Matrix Volume 

^ Total Volume (1-8) 

Again, 0 ^ v^^ £ 1.0. For a material without voids the sum 
of the matrix and filament volume fractions must equal 
unity, 

+ v^ = 1.0. (1-9) 



C. BACKGROUND 

When a tensile test is performed on a composite material 
specimen to determine an ultimate tensile strength the 
value obtained is a macroscopic property. When design 
specifications require certain properties these also are 
macroscopic properties. To accurately determine the behavior 
of a composite and to bridge the gap between design require- 
ments and design practice, macroscopic behavior must be 
related to microscopic behavior within the composite. 

Numerous analytical approaches for the determination of 
composite properties have been utilized. These include 
strength-of-materials , stress function, finite difference, 
and finite element techniques, among others. Chamis and 
Sendeckyj [ 2 ] present a comprehensive review of the various 
analytical methods. 

Specific examples of micro-stress analysis include the 
longitudinal shear loading solutions of Adams and Doner, [s] 
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and Tsai, Adams and Doner [4]. The longitudinal load 
case was considered by Bloom and Wilson [s] using a series 
solution technique. Transverse normal loading was investi- 
gated by Tsai, Adams and Doner [ 4 "] , and Adams and Doner [6] 

using finite difference analyses, and by Foye [?] using a 
finite element method. Lin, Salinas and Ito [s] employed 
finite element techniques to consider the problem of 
combined loading, by superposing the effects of longitudinal, 
transverse, and longitudinal shear loading. Lin et al 
calculated an initial yield surface based on this combined 
loading analysis and later extended their work into the 
plastic region, Adams also considered inelastic 

loading but restricted his study to transverse normal loading 

The combined loading problem has received the least 
attention though it is the most practical from a design stand 
point. It is of great importance to the designer to know 
what combination of macrostresses will result in a micro- 
stress distribution sufficient to cause yield. The initial 
yield surface provides this information. 

D. SCOPE AND OBJECTIVES 

The analytical approach used in this investigation was 
originally developed and programmed by Lin, Salinas, and 
Ito. []8j The author's objectives have been 

i) Modify the program to allow use of an overlay 
system structure in order to reduce the required computer 
storage. 

ii) Substitute the use of random access disk unit for 
the original tape storage methods. 
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iii) Perform a parametric study of the effects of 
constituent properties on the properties of the composite. 

iv) Compare several realistic material combinations to 
evaluate their structural efficiency. 

v) Provide selected initial yield surface plots for 
design assistance. 

vi) Make some general observations relating composite 
behavior to constituent materials. 
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II. METHOD OF ANALYSIS 



The analytical method described in this section was 
developed by Lin, Salinas, and I to [s]. The present discus- 
sion is a general presentation of the basic principles of 
that analysis. Appendix A is a detailed description of 
the method. 

The basic assumptions necessary to the analysis are, 
i) The filament and matrix materials are isotropic 

and homogeneous. 

ii) The filament and matrix materials are perfectly 
bonded . 

It should also be noted that the macroscopic stress 
state is one of plane stress, with components S^, S„, and 

X ^ 

^xz* 

The microscopic stress distribution within a unidirec- 
tionally reinforced composite depends upon the arrangement 
of the filament material within the matrix. In or’der to 
facilitate this analysis the filament array is assumed to 
be a doubly periodic rectangular array. In such an array 
the two-dimensional (x-y) surface may be divided into 
identical segments called basic blocks, of width 2a and 
height 2b as shown in Fig. (2). Away from the edges of the 
composite the microscopic stress distribution should be 
identical within each block. Due to symmetry within a basic 
block only one quadrant need be analyzed, as shown in Fig. (3) . 
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Figure (2). Rectangular Array of Filaments and a Basic Block [sj 



y 




Figure (3) . A Quadrant of the Basic Block [s] 
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In order that no voids be formed and that adjacent 

blocks neither separate or overlap it is necessary that 

the rectangular basic blocks remain rectangular during 

deformation. In order to produce a combined macroscopic 

plane stress state, as well as maintaining the rectangular 

form of the blocks, four specific load cases are defined. 

These four cases, shown in Fig. (4), are superposed to 

give plane stress and boundary compatibility. 

Load cases I and II are plane strain states with 

e =0. Load case III is a generalized plane strain case 
z 

with unit strain. Load cases I, II, and III are superposed 

to give = 0. Load case IV is a state of longitudinal 

shear. All four load cases have = 0. The boundary 

conditions for each problem are shown in Fig. (4). 

These four load cases are solved by the finite element 

method using linear strain triangles. Felippa ^2] presents 

a detailed discussion of finite element theory and a 

derivation of the linear strain triangle. The finite 

element model utilized for the circular filament cross 

section problems is shown in Fig. (5) . This model was 

adapted from that of Lin et al [9]. In this work convergence 

tests were conducted to establish the validity of the model. 

Each load case yields a microstress distribution and 

associated macrostresses. The microstresses are designated 

by • , with subscripts determining the stress component 

I 

and superscripts denoting the load case. For example 
is the microstress at a specific point associated with 
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DofornaLion iltates [B] 




Figure (5) . 



Finite Element Model for Circular Filament 
Cross Section. 
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problem I. Macroscopic stresses are represented by 

with subscripts and superscripts similar to microstresses. 

Appendix A treats the solution of the individual 
problems and demonstrates how the results may be combined 
to relate the combined loading microstresses at a point to 
macrostresses as 



[a] = [A] [S] (ll-l 
sea Appendix A, Eq. (A-13), tor the explicit form of these 
matrices . 

AS noted, the microstresses vary from point to point. 
Therefore the A matrix, evaluated as 



[A] = M [S] 



(II-2) 



is also point dependent. 

in this thesis the Von Mises yield criteria was taken 
to govern initial yield of the composite. In terms of 
microscopic stresses, yielding occurs when [l] 



j = a 

2 xz 



2 + a cr ^ + 



yz 



xy 
3 o 












r^x"] 



= 1 Y ^ 



(II-3) 



where Y is the elastic limit of the material for the tensile 
o 

stress state. 

Using Eq. (II-2) to relate microscopic and macroscopic 
stresses Appendix A shows that the yield equation, Eq. (II 3) 
may be written as 
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fit. 



=1®/ + + C 3 SX + = i (II-4) 

The terms C^, C 2 / C^/ and C^ are second order polynomials 
of the components of the A matrix and are functions of 
position. For given values of S and S , Eg. (II-4) may 
be solved for the value of necessary to cause yield. 

Each point within the fundamental block will yield a dif- 
ferent value of S . Each combination of and S^_ yields 
two values of S_ at each point, a maximum and a minimum. 

The minimum of the maximum values and the maximum of the 
minimum values determine two yield points for the given 
values of S and S . Continuing in this manner the initial 

X X z 

yield surface for a composite subjected to a macroscopic 
plane stress state is calculated. 
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III. EVALUATION OF NUMERICAL RESULTS 



A. INITIAL YIELD SURFACE CHARACTERISTICS 

The two-dimensional plots of the yield loci of this 
investigation have the general form of Fig. (6). Each 
locus represents a specific value of longitudinal shear 
stress, S„„, with the outermost locus corresponding to 
S„„ = 0, as seen in Fig. (6) . 

Each point on a yield locus indicates that some point 
within the material would yield under the given macrostress 
state. The analysis used in this investigation provides 
the location of the area of yield associated with each 
point of a yield locus. These area locations refer to the 
triangular elements of the finite element grid. Fig. (5). 
Figure (6) is the yield surface for a 50% Boron/6061 
Aluminum composite. The elemental areas that first reach 
yield under specific loading conditions are listed in Table 
I. The corresponding macrostress states are noted on Fig. 

(6) . Proceeding in this manner a particular triangular 
element could be associated with every point on each yield 
locus. This information could be useful to a designer in 
determining if failure is associated with a particular 
region. 

On initial consideration it seems somewhat unusual 

that a yield locus is not symmetric about the S s axes. 

X z 

To explain this behavior first recall the yield equation of 
the previous section, Eq. (II-4) . 
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Boron/6061 Aluminum Yield Surface (v 



The exact definition of the terms C 2 / C^, and can be 

ascertained from Eq. (A-22) . 



TABLE I. INITIAL YIELD FOR S = 0 

xz 



Point 


S (ksi) 

X 


S (ksi) 
z 


Yield Location 
(Element) 


A' 


0 


114 


66 


A' 


0 


-114 


66 


B 


+ 15 


+ 112 


67 


B' 


+ 15 


+ 75 


88 


C 


+ 30 


+ 100 


67 


C 


+ 30 


+ .45 


88 



As stated, the points of the initial yield surface are 
calculated by setting to specified values and 

calculating the value of resulting in a microscopic stress 
sufficient to cause yield. Assume that the macrostress 
combination S^, S^j^) is known to be such a yield point. 

It will be shown that for different combinations of S and 
^zl' i*®*/ - S^]^) etc., a specific type of symmetry 

exists . 

For (Sjj, Eq. (III-l) becomes 

+ CjS^iS^ + -I- (III-2) 
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2 



(III-3) 



For (-S^, -S^^) Eq. (III-l) is 

ClSzl“ + C2S^;^S^ + + C^S^2^= I 1^0 

Equations (III-2) and (III-3) are identical and show 
that yield occurs at the same point and for the same 
absolute values of macrostresses in quadrants I and III 
of the yield locus, as shown in Fig. (7). This indicates 
symmetry about the S^„ axis, perpendicular to the (S„,S™) 
plane. 

For combinations of (S„, S ,) such that their product is 

X Z X 

negative, (Sj^, or (~S^, , Eq. (III-l) becomes 

Comparing Eqs. (III-2) and (III-4) it is seen that if 
Sjj, S^l) is a valid solution of Eq. (III-l) as was assumed, 
then (S^z, -Szl) and ”^x' ^zl^ cannot be a solution 

For a positive S and negative S„, solution of Eq. (III-l) 
will yield a different minimum S^, say S^2' This corresponds 
to a different point in the basic block than does the 
solution for quadrants I and III. The combinations (S^., 

Sxz) and (”S^, S^2r identically satisfy Eq. (III-l) 

establishing the equality of yield points in quadrants II 
and IV of Fig. (7). The symmetry of the yield locus about 
the axis is thus established. 

Another interesting phenomenon apparent in some yield 
surfaces is the occurrence of square ends on the yield loci 
corresponding to small values of This is shown in Fig. 
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(8) . Generally this indicates that yield is occurring in 
the filament material rather than the matrix. Note that 
the scale is considerably larger than the scale. 

The longitudinal stress values, S^, in the vicinity of the 
square ends are many times larger than the transverse, S^, 
values. The values are so dominant that increasing the 
transverse stress from point A to point B of Fig. (8) has 
little effect on the yield surface. Once the macrostress 
state B is reached yield occurs in the matrix material. 




As stated the yield locus plot is a two-dimensional 
representation of a three-dimensional surface. Thus the 
separation between successive yield loci is an indication 
of how steep the slope of the yield surface is in that area. 
The dark areas of the plots indicate an extremely steep slope, 
such that adjacent loci actually plot very close to one 
another . 
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(8) . MODMOR I/NARMCO Epoxy (v 








I 






B. EVALUATION OF COMPOSITE MECHANICAL PROPERTIES 

1. Composite Strength 

Composite strength properties were determined 
directly from the yield surface associated with each 
particular composite. 

i) Longitudinal yield strength was taken as 
points S on the initial yield surface when S = Sv= 0. 

ii) Transverse yield strength was taken as points 
on the initial yield surface when S 2 = = 0. 

iii) Longitudinal shear yield strength was taken as 
the point S , on the initial yield surface when S^= S- = 0. 

2 . Composite Elastic Constants 

Composite elastic constants were calculated from the 
results of the individual longitudinal, transverse, and 
longitudinal shear loading problems . The macroscopic 
stresses associated with each loading case were used in 
evaluating the transversely isotropic stress-strain relations 
developed in Appendix B. As discussed in this appendix a 
transversely isotropic material is an orthotropic material 
with a single plane of isotropy. It should be noted that 
the assumption of transverse isotropy limits the applicability 
of this method of determining the elastic constants of a 
composite. Composites are not generally transversely 
isotropic. The particular case treated in this investigation 
is that of a square array of filaments. For this particular 
fiber distribution the composite properties approximate those 
of a transversely isotropic material. For more general 
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filament arrangements, such as a rectangular array, further 
consideration must be given to the applicability of these 
stress-strain relations. The analytical method employed 
to calculate the displacements, stresses and strains is 
not effected by this restriction however. The calculation 
of these quantities is applicable to a general rectangular 
filament array. 

For convenience the transversely isotropic stress-strain 
relations are repeated here. 



a) 


= 


E ^^x 
T 




VL 

El 


b) 


ty 


t ‘"y ■ 

T 


^mS ) 

T X 




c) 


e = 
' z 


- <Sx 


+ Sy) 


+ 1- . 


d) 


^xy 








e) 


Syz 


= — s 
Gl y- 






f) 


e 

xz 


1 

= 7^ S 
Gl xz 







(III-5) 



In the transverse loading case, problem I, the macro- 
scopic strain, corresponding to the unit displacement 

imposed, was calculated. The macroscopic strains Cy and 
e_ are equal to zero. 

In the longitudinal loading case, problem III, the macro- 
scopic strain equal to unity. In this problem 

— ^y 0 • 
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Equations (III-5 ) , evaluated for these two loading 
conditions, were solved for longitudinal and transverse 
composite stiffness and Poisson's ratios, E , E , v , v . 

Li T Ij T 

The longitudinal Poisson's ratio was evaluated directly 
from Eq. (III-5c) in terms of load case I as 






s^+ sj 

X y 



Longitudinal stiffness was evaluated from Eq. 
terms of load case III as 



(III-6) 
(III-Sc) in 



m III III 

E = - V, (S^ + sl^^) 

L ^ 



(III-7) 



L'^'x ■ ^y 

Equation (III-5b) may be evaluated in terms of load cases 
I and III and the resulting equations solved simultaneously 
give 






V, 



si" s^- s”" si 
y s. 



X 



gill. gill gl 



X 



X 



(III-8) 



V 

T 




gi gill- $1 sill 

2 V V 7 

glllgl . gl gill 



X 



X 



(III-9) 



The longitudinal composite shear modulus, G, was obtained 
from the results of the longitudinal shear loading case, 
problem IV. The shear strain corresponding to the imposed 
unit displacement was calculated and used in conjunction 
with the longitudinal shear strain Eq. (III-5f) . 
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3 . Specific Strength and Specific Stiffness 
Specific strength and specific stiffness were 

calculated by dividing the respective properties of the 
composite by the composite density. Densities of the 
various material combinations were calculated from the 
specific weight data of Table II, which lists properties 
for all materials used in this study. 

4 . Stress Concentration Factor 

Three specific values of the stress concentration 
factor are calculated in this investigation, corresponding 
to three specific loading cases. All comply with the 
general definition of the ratio of the maximum microscopic 
stress to the macroscopic stress. 

i) SCFX refers to the transverse loading problem, 
problem I. It is defined as the maximum microscopic stress, 
a„/ divided by the average surface traction, S^. From Eq. 



(1-4) , 



(cfv). 



SCFX = 



max 



(III-IO) 



ii) SCFZ refers to the longitudinal loading problem, 
problem III. It is defined as the maximum microscopic 
stress, a , divided by the average surface traction. From 



Eq. (1-5), 






SCFZ = 



max 



(III-ll) 



iii) SCFXZ refers to the longitudinal shear loading 
problem, problem IV. It is defined as the maximum micro- 
scopic shear stress o„_, divided by S^_, the average shear 

X Z X z 

stress of this problem. 
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SCFXZ 



max 



(III-12) 



S 

xz 

C. COMPARISON OF NUMERICAL RESULTS 

In order to establish the validity of the analysis 
used in this investigation comparisons were made with avail- 
able theoretical and experimental results. It is not the 
intention here to establish the superiority of one method 
of analysis over another but rather to ensure that the 
solutions are compatible. 

1 . Transverse Stiffness, E^ 

Using a finite difference approach Tsai et al [ 4 ]] 
determined the transverse stiffness of a 40% fiber volume 
Boron/6061 Aluminum composite as E,p = 20 x 10® psi. The 
present analysis gives E^ = 20.97 x 10® psi. 

Adams and Doner, [6^ also using a finite difference 
analysis, calculated a theoretical value of transverse 
stiffness for a 47% E-glass/Epoxy model as 1.55 x 10® psi, 
and reported an experimental value of 1.30 x 10® psi. This 
analysis gives 1.42 x 10® psi. Adams and Doner also list 
a comparison between their results for Glass/Epoxy and 
several other analyses, including finite element, complex 
variables, and other solution techniques. These results 
are compared to the present analyses in Table III. 

Transverse stiffness for Boron/Epoxy and Glass/Epoxy 
also fell within the variation bounds calculated by Dow 
et al ( 15 ^ . Further comparisons were made with results of 
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TABLE II. MATERIAL PROPERTIES 





Material ^ 


E (psi) 


V 


Y (psi) 2 


Specific Wt. 
(Ib/in® ) 


1. 


Boron 


60x10® 


0.2 


4.8x10® 


.085 


2. 


S-glass 


12x10® 


0.2 


5.0x10® 


.090 


3. 


E-glass 


10.6x10® 


0.22 


4.0x10® 


. 090 


4. 


MODMOR I 
(graphite) 


60x10® 


0.25 


2.5x10® 


.072 


5. 


MODMOR II 
(graphite) 


40x10® 


0.25 


4.0x10® 


.063 


6. 


THORNEL 25 
(graphite) 


25x10® 


0.25 


1.8x10® 


.052 


7. 


THORNEL 40 
(graphite) 


40x10® 


0.25 


2.5x10® 


. 052 


8. 


6061 Aluminum 


10x10® 


0.3 


.346x10® 


.100 


9. 


2024 Aluminum 


8.1x10® 


0.34 


.130x10® 


.100 


10. 


NARMCO 2387 
(epoxy) 


.5x10® 


0.3 


.300x10® 


. 045 


11. 


4617 Epoxy 


.78x10® 


0.34 


.800x10® 


.045 



1. Materials 1 through 7 are filament materials. 
Materials 8 through 11 are matrix materials. 

2. Y is the elastic limit of the material for the 

o 

tensile stress state. 
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Ekvall |.6] , Whitney ^7] and others. The majority of 
these results are in the form of plotted curves which 
makes direct numerical comparisons difficult but does 
provide a general comparison of results. 



TABLE III. COMPARISON OF TRANSVERSE NORMAL LOADING SOLUTIONS 



Eg/Ejjj=20 Vf=.2 Vjjj=.34 v^=.55 

Analysis Transverse Stiffness, E^ 



Adams and Doner [js] 


1.86x10® 


psi 


Wilson ji;^ 


1.86x10® 


psi 


Pickett (14] 


1.85x10® 


psi 


Foye [7J 


1.83x10® 


psi 


Present analysis (47%) 


1.42x10® 


psi 


Present analysis (60%) 


2.05x10® 


psi 


Present analysis (50% estimate) 


1.81x10® 


psi 



2 . Composite Shear Modulus 

Adams and Doner [s] , using a finite difference method, 
have calculated numerous theoretical values of composite 
shear moduli for specific material combinations. For 
Glass/Epoxy they obtained a shear modulus of .417x10® psi 
and .589x10® psi for 40% and 55% models respectively. The 
present analysis yields .447x10® psi for a 47% model of the 
same material. 
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For Boron/Epoxy Adams and Doner obtained .441x10® psi 
and .658x10® psi for 40% and 50% models, compared to the 
present values of .476x10® psi for a 47% model. 

Direct numerical comparisons were available from 
Blakslee et al who utilized Whitney's micromechanics 

model (l'^ based on an Airy stress function. The particular 
graphite investigated was "Thornel" (a commerical graphite 
produced by Union Carbide) . For Thornel 25 in a 50% model 
Blakslee reported a theoretical value of the composite shear 
modulus of .518x10® psi and an experimental value of .60x10® 
psi. This compares to a present value of .531x10® psi. 

Similar values for a 60% model of Thornel 40 are .931x10® psi 
theoretical and .740x10® psi experimental, compared to a 
present value of .768x10® psi. These comparative values for 
Thornel are summarized in Table IV. 

Graphical comparisons were available with Ekvall ^6] , 
Tsai et al, [ 4 ^ and Whitney p.7] . 

TABLE IV. SHEAR MODULUS COMPARISON FOR THORNEL/EPOXY 





Shear Modulus 


, G(psi X 10"®) 


Analysis 


THORNEL 25* 


THORNEL 40** 


Blakslee (I 8 ] 


.518 


.931 


Experimental [l^ 


.600 


.740 


Present Analysis 


.531 


.768 



* THORNEL 


25 


E£/Em= 


32 


V£“ • 25 


^m=-34 


Vf = 


** THORNEL 


40 


Ef/E = 

t m 


51 


V£~ *25 


V =.34 
m 


^f= 
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The graphical results of Fig. (9) were compiled 
by Reidinger et al jl9] . This compares the values derived 
by Tsai [20] , Dow et al |2^ , and Ekvall |l^ for Glass/Epoxy. 
The results of this analysis are marked by X. Both Tsai and 
Dow et al obtained their results from consideration of 
variational bounds. Tsai however employs a "contiguity" 
factor dependent on the amount of contact between filaments. 
The present analysis and that of Dow et al assume that the 
filament is completely surrounded by matrix material. This 
accounts for the closer correspondence of the present results 
to that of Dow et al. 




E^=10.6xl0® psi Ej^=0.5xlC® psi v^=.22 Vjj^=.35 
Figure (9). Shear Modulus, G, vs. Fiber Volume Fraction,v^ 
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3. Longitudinal Stiffness, 

Longitudinal loading has not been considered as 
frequently as transverse or shear loading, resulting in 
fewer nvimerical comparisons. Blakslee ^s] again gives 
theoretical and experimental values for Thornel/Epoxy 
composites. For the 50% volume Thornel 25, he reports a 
theoretical value of 12.78 x 10® psi, and an experimental 
value of 11.07 x 10® psi. This compares to a value of 
11.96 X 10® psi in this study. For 65% Thornel 40 similar 
results are theoretical 26.20 x 10® psi, experimental 
24.6 X 10® psi and a present value of 24.35 x 10® psi. This 
data is summarized in Table V. 

Reasonable correlation was also obtained in 
comparison with graphical results of Tsai [2(^ and Ekvall [l6] . 

TABLE V. LONGITUDINAL STIFFNESS COMPARISON FOR THORNEL/EPOXY 



Analysis 




Longitudinal Stiffness, 
THORNEL 25* 


E, (psi X 10 ®) 
■^THORNEL 40** 


Blakslee [is] 




12.78 


26.20 


Experimental 


[18] 


11.07 


24.60 


Present Analysis 


12.44 


24.35 



* THORNEL 


25 


E,/E =32 
f m 


Vg= . 25 


^^m=-34 

m 


v^= . 5 


** THORNEL 


40 


Ef/Em=51 


V£= . 25 


Vm=.34 


Vj= . 65 
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4 . Composite Yield Strength 



Strength comparisons were more difficult to obtain. 

In some cases strength quantities described by identical 
terminology and applied to seemingly similar materials dif- 
fered by several hundred per cent. One reason for such 
discrepancies is the lack of consistent definition of 
strength terms. Valid comparisons only have meaning if 
the quantities investigated are precisely defined. 

Another important point is the effect of the matrix 
yield strength on composite yield strength. Matrix mate- 
rials of the same class, such as different brands of epoxy, 
normally have elastic constants that show little variation 
from one another. The difference in yield strength however 
can be considerable. The result of this difference is shown 
in Table VI, comparing a 40% filament model of Boron/2024 
Alumin\am with a matrix yield of 1.30 x 10** psi and Boron/6061 
Aluminum with a matrix yield of 3.46 x 10 psi. 

TABLE VI. EFFECT OF MATRIX YIELD STRENGTH 
ON COMPOSITE PROPERTIES 








Boron/6061 

Al\aminum 


Boron/2024 

Aluminum 


% Difference 


^L 


(psi) 


20.53x106 


18.86x106 


8. 


^T 


(psi) 


31.07x106 


28.71x106 


8. 


G 


(psi) 


6.33x106 


5.28x106 


19. 


S 

L 


(psi) 


94.53x103 


38.9x103 


143. 


"t 


(psi) 


28.50x103 


11.20x103 


154. 
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As seen in Table VI the change in elastic properties 
is insignificant in comparison to the change in strength. 



The change that does occur in the elastic properties is due 
to minor differences in the elastic properties of the two 
matrix materials. The reason- for the very large increase in 
strength may be seen from Fig. (10). 




Figure (10) . Strength Comparison of Different Alloy 
Matrix Materials. 

The slope of curve A represents the composite 



two composites, each having the same filament material 
but with different alloys of the same matrix material. 

As seen in Table VI the values of in such a case are 
quite similar. Curve B represents the stress-strain curve 
of one matrix material, curve C that of the other matrix 
material. They are assumed to have approximately the 
same elastic modulus (as, for instance, two different 
aluminum alloys) but with different yield stress. For 
the case of matrix B the yield stress will be reached at 
a composite strain e^. The resultant composite stress is 



b 




A 



S 



L. 



1 



S 



L. 



Y 



'2 



o 



B 



Y 



o 



c 




S 



L 



1 
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If the matrix material is that corresponding to curve 
C, matrix yield will occur at composite strain resulting 

in composite stress S . Thus a variation in matrix 
strength alone may markedly affect the composite strength. 

No comparison of strength terms should be attempted 
unless the complete properties of both constituents are 
specified. 

The example of Boron/2024 Aluminum of Table VI is 
directly comparable with the results of Adams [li^. For 
identical material properties using a finite element 
analysis with constant strain triangles Adams reported a 
theoretical transverse yield strength of 11,200 psi, 
identical with the results of this analysis. 

Longitudinal strength may. be readily analyzed by 
strength-of-materials methods. Holister and Thomas [ 22 ] 
developed a relationship for determining the ultimate 
tensile strength of a composite undergoing axial tension 
by considering the relation 

= OfVf + a^Vjjj (III-13) 

where v. and v are the filament and matrix volume fractions, 
f m 

respectively. 

For the case of axial load for continuous fibers, as in 
this study, the strain will be constant. Holister and 
Thomas considered the particular case where the strain is 
such that the ultimate stress of the filament is attained 
while the matrix material is still in the elastic range. 
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This is equivalent to assuming a stress-strain relation 
corresponding to Fig. (11) . Curve A represents the 
composite stress-strain relation and curve B that of the 
matrix. 



a 




Figure (11) . Stress-Strain Diagram for Failure Occurring 
in the Filament. 



If the matrix remains in the elastic region, the 
strain fot ultimate strength, must be less than 

that necessary to cause initial yield in the matrix. Under 
these conditions 

^ult = ail-U) 

Equations (III-13) and (III-14) along with the relation 



give 



f m 



= 1 . 



(III-15) 



Suit = Suit + tuiA (III-16) 



Reidinger et al reported an equivalent form of this 
equation as 



ult 



= v.Y 

f o. 



+ V 



m 



'm 



ult 



(III-17) 



49 



where a. 



is defined as the stress in the matrix at the 



ultimate strain of the filament. This is shown in Fig. (11). 

The yield strength of the filament is Yq . 

f 

It must be repeated that these relations are for a 
composite that fails while the matrix is still elastic. It 
should be pointed out however that Holister pursues the 
matter further to consider the case of the matrix being in 
the plastic range. 

For the more common case of the matrix material 

yielding prior to the filament reaching yield, relations 

similar to (III-13) and (III-14) may be utilized to find 

the initial yield strength. In the case of Fig. (12) the 

strain sufficient to cause initial yield for longitudinal 

load, S , is marked by . 

^ o 




Figure (12) . Stress-Strain Diagram for Failure Occurring 
in the Matrix. 

For this case 



Y 




where Y^ is the axial yield of the matrix, 
m 

(III-13) and (III-18) yields 
Sl ' ''m^o 

m 



(III-18) 



Using Eqs. 



(III-19) 
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The stress in the filament at the yield strain of the 
matrix is where 




The results for longitudinal strength calculated in 
this study differed by approximately 2% from those found 
by Eq. (III-20) . 

5. Composite Poisson's Ratio 

a. Longitudinal Poisson's Ratio . The longitudinal 
Poisson's ratio, Vj^, as calculated from the transversely 
isotropic stress-strain gelations, Eq. (III-5) , exhibits a 
linear behavior with respect to filament volume fraction. 
The results of Fig. (13) for E-glass/Epoxy correspond to 
those of a law of mixtures relation. 






(III-21) 



This is in accord with the theoretical results of Rosen et 
al, ( 2 ^ Whitney and Riley, ^?] and Tsai (|o]. 

Specific experimental results were available from 
Blakslee et al [is] for THORNEL/Epoxy . This data also 
showed adherence to a rule of mixtures. 

b. Transverse Poisson's Ratio . The comparison of 
the transverse Poisson's ratio, v^/ with other theoretical 
studies raised several interesting points. Values found in 
the literature ranged from 0.7 to considerably less 

than 0.1. (i^ Thus the literature indicates a range of the 
composite transverse Poisson's ratio that far exceeds that 
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Composite Longitudinal Poisson's Ratio (v 




Figure (13) . 



Longitudinal Poisson's Ratio vs. 



Filament Volume 
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of the Poisson's ratios of the constituents. The higher 

values exceed the theoretical maximum of 0.5 associated 

with isotropic materials. [2^ 

In this investigation using the transversely isotropic 

stress-strain relations, Eqs. (III-5) , results were obtained 

for V . With the addition of a small amount of filament 
T 

material the composite transverse Poisson's ratio takes on 
a value greater than the Poisson's ratio of either constit- 
uent. For the materials considered the increase for a 30% 
filament volume was approximately 22% above that of the 
matrix. In all cases the value of remained below 0.5, 

The transverse Poisson's ratio decreased with increasing 
filament volume. At high filament volume values were 
less than the Poisson's ratio of either constituent. The 
maximum decrease noted in this study was 47% below the 
filament Poisson's ratio. 

To investigate the cause of the transverse composite 
Poisson's ratio having values outside the range of its 
constituents several additional computer analyses were 
made. The constituent and composite properties of these 
analyses are listed in Table VII. 

If the transverse Poisson's ratio was dependent only 
on the Poisson's ratios of the constituents and the filament 
volume, the results of analyses 1, 2, and 3 below would be 
identical. As noted in the previous section, the longi- 
tudinal Poisson's ratio, v , is dependent only on these 
quantities. This fact is borne out by the data of Table VII. 
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TABLE VII. TRANSVERSE POISSON'S RATIO STUDY 



Analysis 




^f 


^f 


V 

m 




^T 


1 


1 


.7 


.25 


.25 


.25' 


.250 


2 


6 


.7 


.25 


.25 


.25 


.219 


3 


120 


.7 


.25 


.25 


.25 


.109 


4 


120 


.7 


.25 


.30 


.26 


.134 



The results of this study indicate that the Poisson's 
ratio in the transverse direction is also greatly dependent 
on the ratio of constituent moduli of elasticity. The 
greater the ratio E^/E^^, the greater the variance from the 
Poisson's ratio of the constituent materials. 

At small filament volume fractions the transverse 

Poisson's ratio exceeds that of either constituent. At 

large values of filament volume content Vrp is considerably 

less than either constituent. At 50% filament volume the 

transverse Poisson's ratio is approximately equal to the 

average of the constituent Poisson's ratios, irrespective 

of the value of E^/E . This is shown in Fig. (14) . 

f m 

Comparison of analyses 3 and 4 of Table VII shows 
that an increase in constituent Poisson's ratio does effect 
v,p but not as significantly as does the ratio E^/E^^. A 20% 
increase in Vj^ resulted in approximately a 20% increase in 
v^. The effect of changing constituents Poisson's ratio is 
limited however by the range of values for available 
materials . 
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^3,4) * Effect of Con. 3 1 i tuent Stiffness on T^rensve^rse 
Poisson's Ratio. 

In the course of attempting to correlate the results 
of this study with other theoretical results several 
alternative analytical methods were considered. Some of 
these were found to give results for Poisson's ratio 
comparable to other extreme values found in the literature. 
These alternative methods were discarded as erroneous. To 
assist in evaluating the various predicted values of v,p 
two of these methods will be discussed here. 

i) If it is assumed that the composite exhibits 
the same stiffness in all directions the stress-strain 
relations become 



e = — (a - v^a ) 
X E X T y" 



£ 



= ^ 1 



(III-22) 



Using the data of the various load cases these equations 
yield values for longitudinal Poisson's ratio, that 

follow a law of mixtures . As noted in the previous section 
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ij'i.'iiimi/' 





this concurs with most theoretical values. The transverse 
Poisson's ratio, however, takes on extremely small values, 
considerably less than .02 for 60% models. This method is 
not correct due to the inaccurate assumption regarding the 
stress-strain relations of the composite. 

ii) It was also attempted to evaluate the transverse 
Poisson's ratio by simulating a tension test. A finite 
element analysis was carried out on the model shown in Fig. 
(15a) . The boundary conditions were 



II 

>1 


0 


on 


y = 


0 


II 

X 


0 


on 


X = 


0 


II 

X 


1 


on 


X = 


a 



This resulted in a displacement along the boundary as shown 
in Fig. (15b) . 





Figure (15) . Simulated Tension Test of a Composite. 

Consequently the Poisson's ratio, v^, varys along the 
boundary. In the problem considered the values ranged 
from approximately 0.15 to 0.7. The error in this method 
is that the rectangular boundaries of the basic block do 
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not remain rectangular, leading to violation of 
compatibility requirements. 

The above result does however suggest a physical 
understanding of the Poisson's ratio found in this 
investigation. Figure (16a) demonstrates how the results 
of this study require the composite to act under a tensile 
load. Figure (16b) shows how the material would react 
if compatible boundary conditions were not enforced. In 
Fig. (16c) the stresses necessary to ensure the required 
compatible boundary conditions are shown. The stresses 
are such as to cause the displacements curve L 2 to 
coincide with The integral of these stresses is zero, 

i.e. , S =0. 

y 

(on boundary) Oy=0 (on boundary) ty j^0(on boundary) 




Figure (16). Transverse Poisson's Ratio. 



6 . Stress Concentration Factor (SCF) 

Before examining the results obtained for stress 
concentration factors some consideration must first be 
given to the value and interpretation of such results. 
The primary value of SCF is in its relationship to the 
strength of a material. For a homogeneous material 
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initial yield will occur at the point of highest micro- 
stress, that is the location of the SCF. For identical 
materials with different SCF, which may occur as a result 
of loading conditions, initial yield will occur at a lower 
stress level in the material with the larger stress con- 
centration factor. 

For composite materials SCF requires further interpre- 
tation. Due to material discontinuity, stress distributions 
will differ in the matrix and filament materials. For this 
reason SCF should be calculated separately for filament and 
matrix regions of a composite. If the matrix is the larger 
no problem is posed, yield will occur in the matrix material. 
If the filament values are larger further investigation is 
necessary to ensure that initial yield is not occurring in 
the filament. If the ratio of filament to matrix stress 
concentration factor is less than the ratio of filament to 
matrix yield strength, yield will still occur in the matrix. 
That is 




m 



< 1 yield in matrix 




simultaneous yield in 
matrix and filament 




> 1 



yield in filament (III-23) 
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Caution should be exercised in interpreting stress 

concentration factors. Assuming that yield occurs in 

the matrix the macroscopic load required to cause yield 

may be related to the SCF as 

Y 



SCF = 



-’m 



(III-24) 



where is the composite yield strength. This may represent 

a longitudinal, transverse, or longitudinal shear yield 

strength. For different materials with identical matrix 

yield strength, Yq , therefore a smaller SCF will require 

m 

a larger macroscopic load, S^, to cause yield. Thus for 
such materials the SCF may be directly related to yield 
strength. An appealing aspect of this result is that for 
composites having identical matrix yield values, SCF offers 
an immediate evaluation of comparative yield strengths for 
the longitudinal, transverse, and longitudinal shear loading 
cases . 

For any composite with initial yield occurring in the 
matrix, regardless of the matrix yield strength, the compos- 
ite yield strength may be calculated for the three loading 
cases as 



S = Y /SCF (III-25) 

m 

As noted above SCF also allows a rapid check on which 
constituent first reaches yield. 

For transverse and longitudinal shear loading the 
maximum SCF usually occurs in the matrix material. For 
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this reason any reference to SCF for these specific cases 
should be assumed to be a matrix property unless specifically 
stated to the contrary. No such generalization is warranted 
for the longitudinal loading problem. 

Several authors have treated stress concentration factors 
extensively. Adams and Doner, [s] ^6] using a finite difference 
technique, have provided a great deal of numerical data for 
both the transverse and longitudinal shear loading problems. 

For the transverse problem with a 55% Glass/Epoxy composite 
Adams and Doner calculated SCFX = 1.86. They also presented 
comparative data from solutions utilizing other analytical 
methods, ranging in value from 1.49 to 1.78. Analyses were 
not performed here for 55% volume composites, however a 47% 
model gave SCFX = 1.49, and a 60% model gave SCFX = 1.70. 

This data is included in Table VIII. 

TABLE VIII. SCF FOR TRANSVERSE NORMAL LOADING 

E_/E =20 v7=T3 V = . 34 Vjr= . 55 

f m f m f 

Method of Solution Stress Concentration Factor, SCFX 



Adams and Doner [6^ 1.86 
Wilson fs] 1,71 
Pickett |4^ 1.49 
Foye [ 7 ] 1.78 
Present analysis (47%) 1.49 
Present analysis (60%) 1.70 
Present analysis (50% estimate) 1.63 
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For the longitudinal shear problem Adams and Doner [s] 
gave results for SCFX of 40% and 55% models as 1.605 and 
1.787 for Glass/Epoxy and 1.649 and 1.854 for Boron/Epoxy. 
As shown in Table IX these values bracket the present 47% 
solutions of 1.66 for Glass/Epoxy and 1.69 for Boron/Epoxy. 

TABLE IX. SCF FOR LONGITUDINAL SHEAR LOADING 



Material 


Filament Volume 


SCFX 


Source 


Glass/Epoxy 


40% 


1.605 


Adams and Doner 


Glass/Epoxy 


47% 


1.66 


Present Analysis 


Glass/Epoxy 


55% 


1.787 


Adams and Doner 


Boron/Epoxy 


40% 


1.649 


Adams and Doner 


Boron/Epoxy 


47% 


1.69 


Present Analysis 


Boron/Epoxy 


55% 


1.854 


Adams and Doner 



61 




I 



IV. PARAMETRIC ANALYSIS OF COMPOSITES 



The purpose of this section is to analyze the effect 
that individual constituent material properties have on 
the mechanical properties of the composite. Properties 
considered include the modulus of elasticity, E, Poisson's 
ratio, V, and yield strength, Y^, of both matrix and 
filament materials. In this study it is assumed that both 
filament and matrix materials are isotropic. 

The goal is not to investigate every possible combination 
of material properties but rather to establish criteria and 
guidance for predicting composite properties associated with 
specific constituent properties. 

The macroscopic properties that are considered as per- 
formance indices include the initial yield strength, the 
specific strength, longitudinal and transverse stiffness, 
shear modulus, stress concentration factors and specific 
moduli. 

To determine how constituent properties effect composite 
behavior, parameter variations were carried out on ficti- 
tious materials. Variations were performed by starting with 
some actual properties , primarily those of a Boron/Aluminum 
or Graphite/Epoxy composite, and independently varying the 
particular property in question. 



62 



A. EFFECT OF FILAMENT VOLUME 



To evaluate the effect of filament volume on composite 
properties several groups of 40%, 50%, and 60% filament 
volume models were investigated. , These included 

1) Boron/NARMCO Epoxy 

2) High modulus graphite (MODMOR D/4617 Epoxy 

3) High strength graphite (MODMOR II)/NARMCO Epoxy 

4) MODMOR 11/4617 Epoxy 

For MODMOR I/NARMCO Epoxy additional models of 30% and 70% 
were also utilized. The material properties for all mate- 
rials are listed in Table II. 

In a composite material subjected to a constant longitu- 
dinal strain the macroscopic stress is related to the micro- 
scopic axial stresses as 

= VfOf + V a_ (IV-1) 

c I £ m m 

where v^ and v represent the volume fraction of filament 
f m 

and matrix materials. An assumption of this study is that 
the filament and matrix are perfectly bonded. For longitu- 
dinal loading the strain will be constant over any transverse 
cross section, thus; 



Cc = Cf = (IV-2) 

Within the elastic region the stress strain relations give 

E^ = v.E. + V E (IV-3) 

L £ £ mm 

The longitudinal stiffness of a composite should be described 
then by a linear relation, the so called law of mixtures. 
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This line passes through the matrix modulus at 0% filament 

volume (E^ = E when = 0) and through the filament 
L in X 

modulus at 100% filament volume (E = E when v^ = 1) , as 

Ij f X 

shown in Fig. (17) for MODMOR I/NARTiCO Epoxy. 

For the case of transverse stiffness the relationship 
is not linear with respect to filament volume. For loading 
in the X or Y direction a constant strain applied at the 
boundary will not result in a uniform strain throughout the 
material. This is due to the change in material composition 
for cross sections perpendicular to the direction of loading. 
At zero filament volume the material has the properties of 
the matrix, that is the composite transverse modulus, E^, 
must equal the matrix modulus, E^^^. At 100% filament volume 
the material will have the properties of the filament and 
the composite transverse modulus, E^, will equal the filament 
modulus, E^. The transverse stiffness, E,p, versus filament 
volume, v^ , curve must go through 

Vj = 0 

and increase monotonically , terminating at 

E^ = Ef at v^ = 1. 

All cases investigated followed this pattern. The results 
for MODMOR I/NARMCO Epoxy are plotted in Fig. (17) . The 
increasing slope of the curve indicates that an increase in 
filament volume has a more significant effect on transverse 
stiffness at higher filament volume percentages. For example, 
for a MODMOR I/NARMCO Epoxy composite a 15% increase in 
filament volume, v^, from 35% to 50% results in a 42% increase 
in composite transverse stiffness, E,p. A similar 15% increase 
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Filament Volume Fraction 

Figure (17). Composite Elastic Properties vs. Filament Volume 
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I f iil 



in from 60% to 75% results in a 75% increase in composite 
transverse stiffness. For composite shear modulus, G, a 
change in cross sectional material properties indicates a 
behavior similar to that of the transverse stiffness. 

In longitudinal loading problem III, the stress con- 
centration factor SCFZ is calculated. The unit strain 
applied results in a uniform stress in the matrix material 
and a uniform stress of a different magnitude in the 
filament material. For the same constituent material and 
identical unit strain these stresses should not vary with 
changes in filament volume since they must satisfy the 
relations 



EfS 


(IV-4) 


^m^ 


(IV-5) 



The resultant surface traction, S , the average value of 

the microstresses and a acting on the cross section, 

^ m 

will change proportionally with the increase in filament 
volume, according to the equation 

S = VfO. + v_a_ (IV- 6) 

2 f f m m ' ' 

Thus as (a^) r the numerator of SCFZ, Eq. (III-ll) , 
max 

remains constant, the denominator increases lineally 
with increasing filament volume, v^ , and the stress con- 
centration factor itself decreases. 

For a homogeneous material the uniform stress throughout 
the entire block would equal the applied stress and SCFZ 
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would equal one. The curve for SCFZ versus filament 
volume should therefore pass through the points 
SCFZ (Vj = 0) =1, and SCFZ (v^ = 1) =1. 

Consider the case of the SCFZ in the filament. The 
addition of a small amount of filament material to a 0% 
model would result in a negligible change in the surface 
traction. That is, for unit strain 

= 1-0 (IV-7) 

The maximum filament stress in this case would be 



= Eft = Ef (IV-9) 

For a small filament volume the filament SCFZ would be 



f 

SCFZ = — = 

cfw, E 

m 



(IV-10) 



Figure (18) for Boron/6061 Aluminum shows this graphically. 
The points at 0% and 100% filament volume are established in 
accordance with the above discussion while the other points 
are results of the analysis. 

For the matrix SCFZ the maximum stress will remain 
constant while the average surface traction increases with 
increase of v^ , as noted in Eq. (IV-6) . For the matrix case 
therefore the curve will originate at SCFZ (v^ =0) =1, 
and decrease to the point where, with little matrix material 
remaining , 

E 

SCFZ = ^ (IV-11) 

^f 
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Recall from section II-B-6 that if the ratio of the 



filament SCF to the matrix SCF is less than the ratio of 

filament to matrix yield strengths, /Y^ , yield will 

f m 

occur in the matrix. This will determine in each individual 
case which of the stress concentration factors discussed 
above is pertinent. 

SCFX and SCFXZ increase directly with filament volume, 
as shown in Fig. (18). These stress concentrations occur 
in the matrix at the interface. 

B. EFFECT OF FILAMENT YOUNG'S MODULUS 

To evaluate the effect of variations in the modulus of 
elasticity of the filament, E^, on the composite mechanical 
properties four computer analyses were made. Each analysis 
used a 50% filament volume model with material properties 
as listed in Table X. 



TABLE X. FILAMENT YOUNG'S MODULUS VARIATION 



Analysis 


^f 

(psixlO ®) 


^f 


5 

°f_ 

(psixlO 


Em 

(psixlO 


^m 


! 

°m 

(psixlO ) 


1 


30 


.2 


00 

• 


• 

O 
1 — 1 


.3 


.3 


2 


40 


.2 


00 

• 


• 

o 
1 — 1 


.3 


.3 


3 


60 


.2 


00 

• 


• 

o 
1 — 1 


.3 


.3 


4 


70 


.2 


• 

00 


• 

o 
1 — 1 


.3 


.3 



68 




Filanent Volume Fraction 

Figure (18). Stress Concentration Factor vs. Filament 
Volume . 
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Figure (19) displays the transverse stiffness, longitu- 
dinal stiffness, and shear modulus of the composite as a 
function of fiber modulus. This figure shows that the most 
significant effect of increasing the filament modulus is 
the linear increase in the composite longitudinal stiffness, 

, This relates with matrix modulus variation and will be 
Xj 

treated in greater detail in the following two sections. 

The relationship between the composite longitudinal 
stiffness, E , and the filament Young's modulus, E^ , is 
linear. In fact for the same reasons as discussed in 
section IV-A for filament volume variation, the relation- 
ship can be accurately described by a rule of mixtures as 

E = v^E. + V E (IV-12) 

L f f m m 

The plot of the composite shear modulus variation of 
Fig. (19) appears to be linear, although it will be shown 
that this variation must be non-linear. Consider that, for 
the combination of materials considered here , v;hen the 
composite is homogeneous E^ = 10 x lO^psi. Hence 

G = G_ = G = 3.84 X 10® psi 
c f m 

The shear modulus curve of Fig. (19) must vary considerably 
from its apparent linear form to pass through this point, 
as shown in Fig. (25). 

As would be expected, considering the linear variation 
of longitudinal stiffness, the longitudinal yield strength 
also increased linearly with filament modulus increase. 
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60 




Filament Young's Modulus, E^CpsixlO"®) 



Figure (19). Composite Elastic Properties vs. 

Filament Modulus of Elasticity. 
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Variation in transverse and longitudinal shear strength 



was negligible for the range of considered, 
C. EFFECT OF MATRIX YOUNG'S MODULUS 



30-75 X 10 



To investigate the effect of variation of the matrix 
modulus of elasticity on the composite mechanical properties 
a series of four analyses were made. All analyses were made 
with a 50% filament volume model with constituent properties 
as shown in Table XI . 



TABLE XI. MATRIX YOUNG'S MODULUS VARIATION 



Analysis 


Ef 

(psixl0“®) 


^f 


(psixlO 


(psixlO“ ® ) 


m 


m 

(psixlO" ® ) 


1 


60. 


.2 


00 

• 


10. 


.3 


.346 


2 


60. 


.2 


• 

00 


7. 


.3 


.346 


3 


60. 


.2 


• 

00 


5. 


.3 


.346 


4 


60. 


.2 


00 

• 


2.5 


.3 


.346 



As established in Eq. (IV-3) the composite longitudinal 
stiffness varys lineally as 



E 



L 



v^E^ + V E 
^ f mm 



Thus the variation of E with matrix modulus is also linear 

Li 

as shown in Fig. (20). From this figure it is seen that the 

variation of E has a greater effect on the transverse 
m 

stiffness and shear modulus than on longitudinal stiffness. 
An increase of 400% in the matrix modulus of elasticity 
results in percentage increases of 



psi . 
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Figure (20). Coiuposite Elastic Properties vs. 

Matrix Modulus of Elasticity. 
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14 % 



ii) = 207% 

iii) G = 222% 

This also resulted in an increase in transverse and longitu- 
dinal shear yield strength. The complementary effect of 
this result with that of the filament modulus variation 
will be discussed in the following section. 

An interesting aspect of the variation of the matrix 
Young's modulus is that, despite the increase in longitudinal 
stiffness, the longitudinal yield strength, Sj^, decreases. 

The reason for this may be seen in Fig. (21). 



cr 




Figure (21). Effect of Matrix Stiffness on Longitudinal 
Yield Strength. 

For two matrix materials with similar initial yield points 
one with the lower Young's modulus will require a larger 
strain to reach yield. Due to the relatively small effect 
of matrix stiffness on composite longitudinal stiffness, E^^, 
assume that E^ is the same for both cases. The larger strain, 
£ 2 / will result in a larger composite yield strength, Sj^ . 
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D. EFFECT OF CONSTITUENT PROPERTY RATIO 



The results noted for filament and matrix stiffness 
variation may be conveniently correlated through use of 
non-dimensional parameters. For example, composite longi- 
tudinal stiffness, E , may be normalized with respect to 

Ij 

filament stiffness, i.e., E^/E^, and plotted against the 
constituent modulus ratio . Figure (22) is such a 

representation. Longitudinal stiffness may be expressed as 






v^rE^ + V^E 

f f mm 



(IV-13) 



Dividing through by E^ and rearranging terms yields 

E^/Ef = Vg (1 - En^/Ef) + E^/E^ (IV-14) 

Equation (IV-14) describes the configuration of the non- 
dimensional longitudinal stiffness curve. The following 
points are significant; 

i) when E^ = E^ , i.e., a homogeneous material, 

E^/E. = E /E^ = 1. (IV-15) 

L f m' f 

ii) when E^^ = 0, i.e., zero matrix stiffness, 

Eb/Ef = v^(l) = v^ (IV-16) 

Equations (IV-15) and (IV-16) were used to determine the 
extreme points of Fig. (22) . The points noted on this 
figure are the results of this investigation for a filament 
volume fraction, v^ = .5. For comparison some other curves 
were drawn using the endpoints determined by Eqs. (IV-15) and 
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(IV-16) . The non-dimensional longitudinal stiffness 
curve for any filament volume must fall in the triangular 
region of Fig. (22) formed by the v^ = 0 curve and the 
V£ = 1.0 curve. 

The comparative effect of filament Young's modulus 
and matrix Young's modulus variations on composite longi- 
tudinal stiffness may be examined by consideration of Eq. 
(IV-13) . Taking partial derivatives of this equation with 
respect to the individual independent quantities yields 

= V (IV-17) 

3E ^ 

m 

and 

3E 

L 

= V. (IV-18) 

3E^ f 

Thus the effect of constituent stiffness variation on longi- 
tudinal stiffness is directly related to the volume fraction 
of the constituent. For the 50% filament volume case this 
is seen to be true as Figs. (19) and (20) show that filament 
and matrix stiffness variation, respectively, result in 
identical slopes of 0.5 when plotted versus longitudinal 
stiffness . 

These comments concerning comparative effects of matrix 
and filament stiffness variation refer to the change per unit 
stiffness of constituent. The actual range of longitudinal 
stiffness due to matrix and filament stiffness variation is 
considerably different however. Consider Fig. (23), a non- 
dimensional longitudinal stiffness curve for some arbitrary 
filament volume fraction, v^. 
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Em/Ef 



Figure (22) . Effect of Constituent Stiffness on 
Longitudinal Stiffness. 
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EL/Ef 




Figure (23) . Non-Dimensional 



E /E. 
m' f 



Longitudinal 



Stiffness . 



For a given constituent moduli ratio, say point A, a 

specific ratio, E /E , results at point A'. At this point 

L f 



El = A'Ef 



If the constituent ratio is increased to point B by increas- 
ing the matrix modulus while holding the filament modulus 
constant then 



E 



B 



E 



B 



B'E^ 

B 



The ratio of longitudinal stiffness for these two hypothetical 
composites is 



E^ B'Ef (E./E.) 

B B' _ ^ ^ 

A’Ej A' (E^/Ej) 



(IV-19) 



Thus variation in the matrix modulus results in a change in 
composite longitudinal stiffness in direct proportion to the 
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In general 



change in non-dimensional stiffness, E /E^ 

the maximum range of non-dimensional stiffness, E /E^, is 

from the value of filament volume fraction, v^, to 1.0. 

The maximum increase in longitudinal stiffness E , due 

L 

solely to a variation of matrix modulus is 

/ \ 




1.0 




max 



For the 50% filament volume case of Fig. (22) this maximum 
increase is a factor of 



/ \ 



B 



\ 



1.0 

.5 



= 2.0 



max 



If the constituent ratio of Fig. (23) is increased to 
point B by decreasing the filament modulus while holding 
the matrix modulus constant, then 



Ef 



B 



< 



E 



f 



A 



E 





A'E. 



D • 



B'E, 



B 



B 



(IV-20) 



Thus the increase in longitudinal stiffness with variation 

of the filament modulus is dependent on the non-dimensional 

stiffness, E /E-, and the ratio of filament moduli, E^ /E^ . 

^ ^ A B 
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For actual filament materials the range for the filament 
modulus is approximately 10 x 10® to 60 x 10® psi. Consider 
a constant matrix material, for example aluminum. For the 
50% filament volume case of Fig. (22) the range of 
corresponding to this range of filament modulus is 50 to 
1.0. For this 50% case the variation in longitudinal stiff- 
ness due to variation in the-^ilament modulus alone is 



Variation of the filament modulus is thus seen to result in 
a much greater range of longitudinal stiffness than does 
the variation of the matrix modulus of elasticity because 
of the greater range of E^ as compared to E^^^. 

Figure (24) is a non-dimensional representation of the 
effect of constituent stiffness, on a normalized 

composite transverse stiffness, E^/E^. Note that the curve 
approaches E^ = 0 as E^/E^ approaches zero. This indicates 
that when the matrix material is extremely soft in comparison 
to the filament that the macroscopic transverse stress 
required to produce a unit displacement of the block boundary, 
as in load case I, is virtually zero. In the limit, when 
E^ = 0, then E^ = 0. Thus it is theoretically possible to 
produce a composite with a value of transverse stiffness 
approaching zero. 

Comparison of Figs. (19) and (20) shows that the slope 
of the E vs. matrix stiffness curve is greater than that 
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Figure (24) . Effect of Constituent Stiffness on 
Transverse Stiffness. 
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of the E vs. filament stiffness curve. Therefore on a 
T 

unit stiffness basis, variation of matrix Young's modulus 
is more effective than variation of filament Young's modulus 
on the composite transverse stiffness. 

By considering a non-dimensional transverse stiffness 
curve and following the same technique as for the non- 
dimensional longitudinal stiffness the range of the trans- 
verse stiffness due to matrix Young's modulus variation, 
analogous to Eq. (IV-19) , is 



/ \ 










(E^/E ) 


B 




B 






(E^/Ef ) 


max 


A 



max 



(IV-21) 



The range of E^/E^ , as is seen in Fig. (24), is zero to one. 
Thus the theoretical maximum increase of transverse stiff- 
ness due to matrix Young's increase is 



E, 



Er 



■B 



1.0 

0 



^/max 

The maximum variation in transverse modulus due to variation 
of the filament stiffness alone is 



Eip 

A 




(E /E ) E 

A 


Em 








max 


B B 



max 



(IV-22) 



For common filament materials the range of filament modulus 
is approximately 10 x 10® to 60 x 10® psi. Consider a 
constant matrix material, for example aluminum. For the 50% 
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filament volume case of Fig. (24) the range of E^/E^ 
corresponding to this range of filament modulus is .349 
to 1.0. Thus the maximum variation in transverse stiff- 
ness due to variation of the filament modulus alone is 

/ \ 




max 




2.094 



The range of transverse stiffness due to matrix Young's 

modulus variation is greater therefore than that due to 

filament Young's modulus variation. 

Figure (25) is a plot of the normalized composite shear 

modulus, G/Gj, versus the ratio of constituent shear moduli, 

G /G^. This curve substantiates earlier contentions that 
m' f 

shear modulus variation with respect to constituent moduli 
is non-linear. Again this curve approaches zero as the 
constituent moduli ratio approaches zero. 

Again comparing Figs. (19) and (20) it is seen that 
the slope of the composite shear modulus versus matrix 
stiffness is greater than the slope of the composite shear 
modulus versus filament stiffness curve. This indicates 
that, per unit of constituent stiffness, the matrix material 
has a greater effect on composite shear modulus. 

By considering a non-dimensional shear modulus curve 
and following the same technique as for the non-dimensional 

longitudinal stiffness, the range of shear modulus variation 
due to matrix shear modulus variation, analogous to Eq. (IV-19) , 
is 
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(IV-23) 



\ 




G„ 




(G/G.) 


B 




B 


G 




(G/G ) 


^1 


max 


f A 



max 



The range of G/G^ , as seen in Fig. (25), is from zero to 
1.0. Thus the theoretical maximum value of Eq. (IV-23) is 









_ 1.0 



max 



If the ratio is increased by changing the filament 

modulus the maximum range attainable is 



\ 










(G/Gi) Gi 

A A 


S 


max 


(G/Gi) G^ 






B B 



(IV-24) 



max 



For common filament materials the range of shear modulus is 
approximately 4.0 x 10® to 25.0 x 10® psi. For a constant 
matrix material, for example aluminum, the range of G/G^ of 
Fig. (25) corresponding to this range of filament shear 
modulus is .31 to 1.0. These values are for a 50% filament 
volume case. For this case the Eq. (IV-24) becomes 



_A 



max 



1.0 



\ 




1.94 



Variation of the matrix shear modulus thus results in 
a greater range of composite shear modulus than does 
variation of filament shear modulus. 

The effect of Poisson's ratio will be discussed in the 
next section. 
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Figure (25) . Effect of Constituent Shear Modulus on 
Composite Shear Modulus . 
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The primary advantage of the non-dimensional represen- 
tations is that valid comparisons may be made between 
composites without necessarily holding the Young's modulus 
of one constituent constant. 

It should be noted that the non-dimensional plots 
presented here are for a 50% filament volume model. Dif- 
ferent filament volume model's will result in different 
curves . 

E. EFFECT OF POISSON'S RATIO 

Two series of computer analyses were made for evaluation 
of the effects on composite behavior due to change of Poisson's 
ratio of the filament, v^. Four analyses were made using a 
fictitious material approximating the properties of Boron/ 
Aluminum while varying the filament Poisson's ratio from .15 
to .30 in increments of .05. 

Another series of analyses was conducted with MODMOR I 
graphite and epoxy. This was done because Poisson's ratio 
for graphite filaments is not accurately specified and it 
was desired to determine whether or not this would signif- 
icantly effect the results of this investigation. 

In both cases no significant effect was noted on any 
composite mechanical property. All observed variations 
were less than 1%. 

The effect of the variation of matrix Poisson's ratio 
was evaluated through a series of three analyses. Constituent 
properties similar to Boron/6061 Aluminum were used while 
varying the matrix Poisson's ratio from to .30 to .35. Again 
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no significant effect on composite mechanical properties 
was noted. It is expected however that changes in Poisson's 
ratio will produce more significant changes for post yield 
behavior. 

F. EFFECT OF MATRIX YIELD STRENGTH 

To investigate the effect of the matrix yield strength, 

Yq , on composite mechanical properties, a series of four 
m 

computer analyses were made. All used a 50% Boron filament 
volume model with constituent properties as listed in Table 
XII. 



TABLE XII. MATRIX YIELD STRENGTH VARIATION 



Analysis 


^f 

(psixlO"®) 


^f 


Yof 

(psixlO"®) 


E 

m 

(psixl0“®) 


^m 


Yo 

m 

(psixl'O"®) 


1 


60. 


.2 


4.8 


10. 


.3 


.346 


2 


60. 


.2 


4.8 


10. 


.3 


.230 


3 


60. 


.2 


4.8 


10. 


.3 


.130 


4 


60. 


.2 


4.8 


10. 


.3 


.030 



The yield point of a material has no effect on the 
elastic properties of that material. It merely defines the 
limit of elastic behavior. Thus the variation of matrix 
yield strength had no effect on the composite stiffness or 
shear modulus. This is as expected considering the compar- 
ison in section III-B-4 of two aluminum alloy matrix 
materials . 
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Figure (26) shows hypothetical stress-strain curves for 

the composite (curve A) and the matrix (curve B) . For a 

given matrix yield strength, say Yq , a specific composite 

m 

macroscopic yield strength, S^, is defined. If, as stated, 

the yield point of a material does not effect the elastic 

properties, then the stress-strain curves remain the same, 

independent of 'the yield strength variation. For a given 

increase in matrix yield strength, say a 50% increase to 

Yq , plane geometry shows that an identical increase, on 
m 

a percentage basis, must occur in the composite yield 
strength. Figure (27) confirms this hypothesis. 




Figure (26) . Effect of Matrix Yield Strength Variation. 

It should be noted that the above discussion assumes 
that yield occurs in the matrix, as is true for most, but 
not all, composites. 

The implication of this phenomenon is that there is no 
variation in the shape of the initial yield surface due 
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Longitudinal Yield Strength, S (psixlO 




Matrix Yield Strength, Yq (psixlO 

ra 



Figure (27). Effect of Matrix Yield Strength on 
Composite Yield Strength. 
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solely to variation in the matrix yield strength. That is, 
the initial yield surface changes in size only. This can 
be seen in Figs. (C-1) and (C-2) of Appendix C. 



G. EFFECT OF FILAMENT CROSS SECTION 

The majority of the computations included in this investi- 
gation are based on longitudinal filaments with a circular 
cross section. This was done because it is assximed that this 
particular cross section would be the easiest to produce. The 
question arises as to the effect change of shape of filament 
cross section would have on composite properties. Such 
changes might be the result of fabrication. 

To provide an answer to this question two elliptical 
cross sectional filaments were investigated. These were 
identical ellipses with major/minor axes ratio of 1.25 to 
1.0. One model had the major axis oriented along the x-axis 
and one along the y-axis, as seen in Fig. (28) . 



The rule of mixtures and the comments of section IV-A 
leading to its use are independent of cross section shape. 
Elliptical and circular cross section filaments of the same 
volume should therefore produce identical results for the 





X axis Major Axis y a 

Figure (28). Elliptical Cross Sections. 



y axis 
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longitudinal composite stiffness, E . In the three cases 

Xi 

investigated, 50% models of Boron/6061 Aluminum, Boron/NARMCO 

Epoxy, and S-glass/NARMCO Epoxy, this was found to be true. 

The variation in E was less than 1% in all cases, as may be 

L 

seen from Table XIII. 



- TABLE XllX. EFFECT OF -FILAJIENT CROSS SECTION 
ON LONGITUDINAL STRENGTH AND STIFFNESS 



Circular Elliptical 

E (psixlO’®) S (psixlO”'*) E (psixld”®) S (psixldT**) 



Boron/6061 Aluminum 


37.08 


11.4 


37.25 


10.1 


Boron/NARMCO Epoxy 


28.23 


16.71 


28.57 


14.5 


S-glass/NARMCO Epoxy 


6.05 


3.47 


6.105 


3.04 



The transverse composite stiffness reacted as expected. 

In all cases E^ increased in the direction of the major el- 
liptical axis and decreased in the direction of the minor 
axis as noted in Table XIV. This seems reasonable since more 
of the stiff filament material is oriented in the direction 
of the major axis. On a percentage basis the increase in 
transverse stiffness in the direction of the major axis was 
approximately twice that of the decrease in the direction of 
the minor axis. 

An interesting result of the cross section variation 
that is evident in Table XIV is that the transverse strength 
changed in a manner opposite to that of the transverse stiff- 
ness. That is, the transverse strength increased in the 
direction of the minor axis and decreased in the direction 
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of the major axis. This result may be understood by 
considering the micromechanics of the problems in question. 
For ease of discussion call the transverse load case with 
the major axis oriented along the x-axis load case lA, and 
that with the major axis along the y-axis as load case IB. 
These load cases are shown in Fig. (29) . Using this 
notation, the -analysis resul'ts require, for transverse 
normal loading, at yield 



S 





(IV-25) 



The definition of macrostresses requires that the integral 



of the microstresses along the boundary be related as 




(IV-26) 



TABLE XIV. EFFECT OF FILAMENT CROSS SECTION 
ON TRANSVERSE STRENGTH AND STIFFNESS 



Circular 


Major 

^T 


Elliptical 

Axis Minor Axis 

^T ^T ^T 


Boron/6061 Aluminum 


24.4 


3.0 


26.6 


3.5 


23.0 


2.5 


Boron/NARMCO Epoxy 


1.64 


,24 


2.05 


.19 


1.46 


.30 


S-glass/NARI»lCO Epoxy 


1.52 


.245 


1.83 


.195 


1.38 


.295 


E^ (psixlO'®) 


S<p (psixlO""*) 











Figure (29) shows the stress and strain distribution on 
the x-boundary (x = a) for a unit displacement for load 
cases lA and IB. In both cases yield occurs in the matrix 
at a point near the interface, along the x-axis. Both cases 
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Figure (29) . Stress and Strain Distribution on the Boundary of a Transversely 
Loaded Elliptical Filament. 



also show little variation in stress or strain along the 
X-axis between the interface and the boundary. For example, 
for load case A, i.e., major axis along the x-axis, problem 
I for a 50% S-glass/Epoxy model gives a stress, Oy^, at the 
interface of .664 x 10® psi. At the boundary (x = a) 

equals .688 x 10® psi. The same material for load case B 

gave an interf-ace value of cr^ as .366 x 10® psi and a 

boundary value of .389 x 10® psi. Thus the strain to cause 

yield will be approximately the same at the yield point, i.e., 
the interface, and at the boundary (x = a, y = 0) . If the 
strain and stress curves of Fig. (29) are scaled to give 
identical strain values at the boundary (x = a,y = 0) the 
more uniform distribution of case IB will result in a larger 
value of the integral of the stress. This confirms Eqs. 
(IV-25) and(IV-26). 

It seems therefore that the stiffness of the composite 
depends on the degree of filament material alignment in a 
particular direction while the strength depends on the 
maximum filament cross section in a direction perpendicular 
to that of the load. The inference of this is that the 
optimum cross-section shape would approximate that of a 
plane whose orientation depends on the desired properties. 

This is shown in Fig. (30). The filament shape can not be 
a rectangle as the square ends would induce prohibitive 
stress concentrations. 
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Figure (30) . Approximate Optimum Cross Section for 
Particular Properties. 

In the case of circular cross section filaments the 
transverse properties are the same in the y-direction as 
in the x-direction. Thus transverse yield strength in 
the y-direction would be numerically equivalent to that in 
the x-direction. This is shown in Fig. (31a) . Again assume 
the convenient designation of load case A as that of the 
elliptical cross section with the major axis along the 
x-axis, and load case B with the major axis along the y- 
axis. The results of these two load cases may be interpreted 
as transverse strength in different directions of the same 
model as shown in Fig. (31b) . 

On a percentage basis the change in transverse strength 
due to cross section variation was considerable, ranging from 
15% for Boron/Aluminxam to 23% for Boron/NARMCO Epoxy. Of 
more significance however is the fact that the vector sum 
of and Sy is independent of the filament cross section. 
This indicates, in vector notation, for Fig. (31) 

S = S + S = + s' (IV-27) 

c ^c e ^e 
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Or, alternatively. 



is. 



+ S 

Y 



2 \ _ 



(S^ " + s 

e ^e 



(IV-28) 





'HaK' 





= S 



X, 



X, 



^ s. 



Elliptical Filament, Identical Vp 
(b) 

Figure (31) . Transverse Strength Relation to Filament 
Cro5j;s Section. 

This indicates that any increase in strength in one direction 
is directly offset by a decrease of the strength in the per- 
pendicular direction. 

In each case there was a noticeable increase in the stress 
concentration factors SCFX and SCFXZ , although SCFZ remained 
relatively constant. The net result was a reduction in the 
longitudinal yield strength of the elliptical cross section 
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compared to the circular cross section. The decrease was 
approximately 12% in each case. The implication is that a 
circular cross section is the most effective, with respect 
to a given filament volume, in terms of longitudinal strength. 
Adams and Doner [s] reported that stress concentration factors 
are directly related to the maximum radius of curvature of 
the filament. -For a given filament volume a circular cross 
section results in the smallest possible maximum radius of 
curvature. This also implies that the circular cross section 
is the most efficient. 

Another significant result evident in the analysis of 
the elliptical filament is the effect of filament shape on 
the composite shear modulus, G. Normalized with respect to 
the circular composite shear modulus the results obtained 
are listed in Table XV. 



TABLE XV. EFFECT OF CROSS SECTION VARIATION 
ON COMPOSITE SHEAR MODULUS 



Composite 


G, 

Major 


sll/^c ^f 

Axis Orientation 




X 


y 


Boron/6061 Aluminum 


1.345 


.555 


Boron/NARMCO Epoxy 


1.21 


.326 


S-glass/NARMCO Epoxy 


1.15 


.398 



The shape of the filament cross section could be used 
to good advantage by a designer with detailed knowledge of 
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the expected loading of the structure under design. Such 
large variations in some composite properties due to 
relatively small departures from a circular cross section 
also make it extremely dangerous to load a "circular" 
filament near its expected transverse or shear strength. 
Small manufacturing errors can lead to large strength 
discrepancies .r- 
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V. COMPARISON OF MATERIAL COMBINATIONS 



The preceeding sections have investigated the effect 
of various constituent properties on the properties of 
the composite. The goal was to be able to predict the 
comparative properties of various actual material combina- 
tions. In this section several functional material com- 
binations are analyzed. Theoretical mechanical properties 
are compared and the results are correlated with the guide 
lines developed in the preceeding sections . 

Forty, fifty, and sixty percent filament volume models 
are compared for the following combinations of materials, 

i) MODMOR II/NARMCO Epoxy 

ii) MODMOR 11/4617 Epoxy 

iii) MODMOR 1/4617 Epoxy 

iv) Boron/6061 Aluminum 

For MODMOR I/NARI4C0 Epoxy and E-glass/Epoxy , 30% and 70% 
filament volume models were also considered. Single models 
of S-glass/NARJlCO Epoxy, THORNEL 25/4617 Epoxy, THORNEL 
40/4617 Epoxy, and Boron/2024 Aluminum were also investigated 
Primary measures of structural efficiency of a material 
are the specific strength and specific moduli of the material 
These are defined in section I-B, however their importance 
necessitate further explanation. 

In general, design of a particular structure requires 
certain material properties of strength and/or stiffness. 

In modern technology the weight of a structure has gained 
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equal importance with mechanical properties. It is 
significant to note that all common structural materials, 
when considered on a stiffness to weight basis, are remark- 
ably similar. Wood, steel, aluminum, and titanium fall 
within the range of 92.0 x 10 ® psi/(lb/in^) to 
105.5 X 10® psi/(lb/inM for specific Young's modulus. 
Similarly for specific ultimate strength these common 
materials range from 4.60 x 10® to 8.33 x 10® psi/ (Ib/in® ) . 
This data is summarized in Table XVI . Thus for a composite 
to be valuable from a structural standpoint it must exceed 
these values for specific mechanical properties. 

TABLE XVI. SPECIFIC ULTIMATE STRENGTH AND SPECIFIC 
YOUNG'S MODULUS FOR COMMON STRUCTURAL MATERIALS 



Material 


Specific Ultimate 
psi/ (Ib/in® ) xlO 


Strength Specific Modulus 
® psi/(lb/in®)xl0“® 


Steel 


5.16 


105.5 


Aluminum 


6.64 


104.2 


Titanium 


8.33 


97.2 


Magnesium 


5.00 


101.6 


Wood 


4.60 


92.0 



Figure (32) shows the transverse stiffness, E^, and 
the specific transverse stiffness, E^/p, where p is the 
composite specific weight. These results are for the 50% 
filament volume models considered. In the preceeding 
sections transverse stiffness was found to depend primarily 
on the matrix modulus of elasticity. The results of Fig. 
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(32a) show that the transverse stiffness of materials with 
similar matrix properties showed little variation. For 
example, those composites with a matrix of NARMCO epoxy 
had comparable values of transverse stiffness regardless 
of filament properties. The aluminum matrix was so much 
stiffer than any other matrix material studied that the 
performance of" Boron/Aluminum far exceeds that of the 
other composites. 

Figure (32b) shows the specific transverse stiffness. 

Note that the smallest value of composite specific transverse 
stiffness is approximately 80% greater than that normally 
associated with structural materials. 

Figure (33) shows the composite transverse strength, S^, 
and specific transverse strength, S^/p . The constituent 
property that showed the most influence on transverse 
strength was the matrix modulus of elasticity. Figure (33a) 
reinforces this judgement by showing negligible variation 
in transverse strength among those composites utilizing 
identical matrix materials. Progressively stiffer matrix 
materials produced greater transverse yield strength regard- 
less of the filament material. 

Figure (33b) shows the specific transverse yield 
strength for the composites considered. Note that not 
all of the material combinations exceed the stated value 
of 8.33 X 10 ^in. for common structural materials. One 
reason for this is that the quoted values for common 
structural materials are ultimate strengths rather than 
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Figure (32) . Transverse Stiffness Comparison. 
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Figure (33) . Transverse Strength Comparison. 
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yield strengths. Another reason is that the composite 
lamina are designed such that their primary strength is 
in the longitudinal direction whereas most standard struc- 
tural materials are considered isotropic. 

Figure (34) is a plot of longitudinal stiffness and 
specific longitudinal stiffness. It has been shown in 
this study' that composite longitudinal stiffness can be 
accurately predicted by a law of mixtures. For the 50 o 
case of Fig. (34a) the composites should be in order of 
increasing sum of matrix and filament moduli of elasticity. 
This is seen to be true. 

Figure (34b) shows the specific longitudinal stiffness. 
Note that the graphite/epoxy combinations are considerably 
more efficient than the Boron/Aluminum combinations despite 
the large stiffness advantage of aluminum over epoxy. 

Figure (35) shows the longitudinal yield strength and 
specific longitudinal yield strength of the various 
composites . 

The relative longitudinal yield strength of a composite 
is somewhat more difficult to predict than are the other 
properties. Increasing the filament modulus or the matrix 
yield strength has been shown to increase the composite 
longitudinal yield strength. Increasing the matrix modulus 
of elasticity decreases longitudinal yield strength. There 
is no set method however for predicting the combined effect 

of such changes. 
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Figure (34) . Longitudinal Stiffness Comparison. 
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Figure (35) . Longitudinal Strength Comparison. 
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The MODMOR II (graphite) /Epoxy models indicate a 

case where the matrix yield strength was of greater 

significarj.ce than the matrix modulus increase. The 4617 

Epoxy had a greater modulus, which should decrease the 

yield strength, and a greater matrix yield strength which 

should increase the composite longitudinal strength. The 

net result was- a marked increase in composite yield 

strength. Considering the variations individually, as 

in Fig. (36) , this indicates that the increase in matrix 

yield strength, say from ‘ the composite 

j m 2 m 

yield strength from to Sj^. This is shown in Fig. (36) . 

The decrease in matrix modulus results in a composite stress- 

strain curve with a smaller slope. For the same matrix 

yield strength, Y the composite yield is S^. The results 

m 

of the analysis shows that 



S 



1 

L 





(V-1) 



a a 





Figure (36). Combined Effect of Matrix Modulus and 
Matrix Yield (M.ODMOR II/Epoxy) 



107 



1 




The two boron analyses represent the opposite case. 

The change in matrix modulus prevails over the change in 
matrix yield strength as the epoxy composite has a yield 
strength considerably larger than that of the Boron/Aluminum 
composite . 

It should be noted that the longitudinal yield strength 
of the graphite composite MO'DMOR 1/4617 Epoxy is considerably 
less than that of the MODMOR 11/4617 Epoxy. This is signif- 
icant in that the Young's modulus of MODMOR I is 50% greater 
than that of MODMOR II. This apparently contradicts the 
guidelines presented in this study. The reason for this is 
that the conclusions of this investigation are based on the 
assumption that initial yield occurs in the matrix rather 
than the filament. Due to the relatively low yield strength 
of MODMOR I initial yield actually occurs in the filament. 
This is predicted by considering the stress concentration 
factors as outlined in section III-B-6. Explicitly, for 
MODMOR 1/4617 Epoxy, 



Filament yield strength 2.5 x 10® psi 
Matrix yield strength .08 xlO® psi 



Filament SCFZ 



2.106 



Matrix SCFZ 



0524 




(V-2) 



f 2.5 




(V-3) 



m 



108 









i 



(V-4) 



SCFZ^ Y„ 
if o 



SCFZ 

m 

Which meets the criteria of section III-B-6 for initial 
yield occurring in the filament. 

Figure (35b) represents the specific longitudinal yield 
strength. The_ unshaded area, of the bar chart represents 
the specific yield strength from the yield calculations of 
this analysis. The shaded region of the figure represents 
an estimate of specific ultimate strength. The ultimate 
strength was estimated by use of the equation 



ult f Oj m 



As discussed in section III-B-4 o 



m 



(V-5) 

is the stress in the 



ult 

matrix at the ultimate strain of the filament. Again it 
should be noted that this relation assumes that the matrix 
remains elastic up to the failure of the filament. Generally 
this is not true. Using the method of Lin, Salinas, and 
Ito [ 9 ^ the theoretical ultimate tensile strength of a 50% 
Boron/Aluminum composite was found to be 243,000 psi. This 
compares with an estimate using Eq. (V-5) of 305,000 psi, a 
difference of greater than 20%. While the discrepancy is 
rather large the results are considered to be worthy of 
inclusion as a relative measure of the range between yield 
and failure. 



For MODMOR I, MODMOR 11/4617 Epoxy and THORNEL 25 (all 
graphite filaments) yield actually occurred in the filament. 
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This is the reason that the specific ultimate and specific 
yield strength are virtually coincident. 

Also of importance in a manufacturing process is the 
cost of the materials. Table XVIII shows strength per 
unit cost data for selected materials. Cost data was 
obtained from Refs. [2^ and [26] . The extremely high 
strength to cost ratio for g^lass filament serves to explain 
why glass composites are in such wide use despite their 
relatively low transverse yield strength. The present 
interest in graphite filaments can certainly be justified 
by the strength/cost increase anticipated in the near future. 



TABLE XVII. COMPARATIVE STRENGTH PER UNIT COST 



Filament 


Epoxy Matrix v^ 

Strength/Unit Cost 


= 0.5 

(in-lb/dollarxlO""*) 


Boron 


.98 




Graphite 


2.01 




Graphite (1974) 


4.02 




Glass 


42.48 





Numerous other factors also enter into a material 
selection process. Ease of handling and fabrication, 
special processing requirements, and specific desirable 
properties such as thermal expansion may be such factors. 
Consideration of these factors is beyond the scope of this 
investigation . 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

In this investigation a finite element analysis of the 
microscopic and macroscopic stresses of a composite material 
has been carried out. The stresses were used to calculate 
the initial yield surface for unidirectional composite 
materials. The results of these analyses were then used 
to determine the mechanical properties of a laminate of 
composite material. Parametric studies of the effect of 
various constituent properties on the composite mechanical 
properties were then carried out. Table XVIII summarizes 
the results of these parameter variations. 

These results are based on the assumption that initial 
yield takes place in the matrix material. The effects noted 
in this table are those associated with a change of the 
specified constituent property only, that is, with all 
other constituent properties held constant. Caution must 
be exercised in attempting to predict the effects of changes 
in more than one constituent property. 

Other results included; 

i) Composite Transverse Poisson's Ratio, v^. This 
value exceeded that of either constituent. By adjusting 
the ratio of filament Young's modulus to matrix Young's 
modulus, E^/E^, the designer may achieve extreme values of 
transverse Poisson's ratio. 
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TABLE XVIII. EFFECT OF CONSTITUENT PROPERTIES ON COMPOSITE PROPERTIES 
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ii) Longitudinal Initial Yield Strength, Sj^. For 
initial yield occurring in the matrix material the longi- 
tudinal yield strength may be. predicted by 



S = VO, + 

m 

where Y is the yield strength of the matrix material, 
m 

is the stress in the filament at the yield strain of 

the matrix, and and are filament and matrix volume 

I m 

fractions, respectively. 

iii) Initial Yield Location. The computer program 
utilized in this investigation gives results which include 
the location of the yield point for any macroscopic stress 
state. For load case III it may be determined whether the 
composite fails in the matrix or the filament by consider- 
ing the stress concentration factor in the longitudinal 
direction, SCFZ. If the ratio of SCFZ in the filament to 
that in the matrix is less than the ratio of filament yield 
strength to matrix yield strength then yield will occur in 
the matrix. That is, if 



SCFZf/SCFZ < Yq /Yq 

f m 

yield will occur in the matrix. This determines the appli- 
cability of the longitudinal yield equation noted above and 
does not require the application of any particular yield 
criteria. 

iv) Filament Cross Section Effects. Variations from a 
circular filament cross section have significant influence 
on composite properties. Shear modulus and transverse 



113 








I 




properties change considerably. Due to the distribution 
of microscopic stresses on the block, boundaries an el- 
liptical cross section has greater transverse strength in 
the direction of the minor axis. 

B . RECOMMENDATIONS 

i) The logical extension of this investigation is 
consideration of inelastic loading. The analysis used 
in this study is a special case of the elastic-plastic 
analysis of Lin, Salinas, and I to. ^9] 

ii) Further effort seems warranted in the evaluation 
of filament cross section effects. Through use of cross 
sectional shapes and filament distribution it should be 
possible to exert control over transverse and shear 
properties . 

iii) The wide variation in the literature of theoretical 
transverse Poisson's ratio certainly suggests an urgent 
need for accurate experimental results as well as a 
comprehensive survey of analytical techniques. 
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APPENDIX A: YIELD LOCUS CALCULATIONS 



The analytical method employed in this investigation 8 
requires that the macroscopic stress state of the composite 
be one of plane stress. It is also required that the rec- 
tangular form of the basic block. Figs. (2) and (3) , be 
maintained under deformation. To ensure these results 
four distinct loading cases are superposed. These load 
cases and their boundary conditions are shown in Fig. (4) . 

Load cases I and II are plane strain states with e =0. 

Load case III is a generalized plane strain case with 
e = 1.0. Load cases I, II, and III, are superposed to 
give Sy = 0 . Load case IV is a state of longitudinal shear. 
All four cases have = 0. Rectangular shapes of 

the basic block result from all load cases. 

Using finite element methods with linear strain triangles 
1 2] the four loading cases are individually analyzed to 
obtain the microstresses and their associated macrostresses. 
The microstresses are designated with a subscript denoting 

the stress component and a superscript for the load case. 

For example, for load case I the microstresses obtained are 



„I „I -I 
®xy' ^z 



(A-1) 

Similar results are obtained for problems II and III. The 



macrostresses for load case I are designated 
and similarly for load cases II and III, 



(A- 2) 
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Load case IV yields microstress results for 



xz 

and macrostresses 



and 



alV 

yz 



(A-3) 



and (A-4) 

In the notation of Lin, -Salinas and I to [s] the micro- 
stresses of a combined loading state may be related to the 
macrostresses as 



and 



{a} = [a] {S} 
4x1 4x3 3x1 



{a} 

2x1 



[a] {S} 

2x2 2x1 



Where 

T 

= < o a a o > 

X y xy z 

{S =<S^SyS^ > 
_ 

{a} =< 

{S } = < Syz > 



[^J 



^1 


^12 


^13 


^21 


^22 


^23 


^31 


^32 


^3 


A41 


^42 


^3 



(A-5) 

(A-6) 

(A-7) 
(A- 8) 

(A- 9) 
(A-10) 



(A-11) 
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The matrix A may be related to the results of load 
cases 1/ II and III as 



I II 


III ~ 


- 








^x ^x ‘ 


^x 




All 


Ai2 


A13 


I II 


III 




A 


A 


A„,, 


Oy Oy 


"y 


- 


21 


22 


23 


I II 


III 




A 


A „ 


A„, 


*^xy *^xy 


*^xy 




31 


32 


33 


I II 


III 




A 


A 


A,„ 






1 


41 


42 


43 






Sy Sy Sy 






(A-13) 



(A-14) 



For convenience this may be written 

[t] = [a] 

4x3 4x3 3x3 

The matrix A can then be evaluated as 

[a] = [t] M - ■ 

4x3 4x3 3x3 

Thus the combined microstresses, [o] / due to any unifo 
transverse or longitudinal deformation may be evaluated by 
substituting the results of Eq. (A-15) into Eq. (A-5) as 



(A-15) 



(A-16) 

In a similar manner for load state IV and an analogous 

Load case V associated \ 

A may be evaluated as 



{a} ^ [t] [C' 

4x1 4x3 3x3 3x1 



load case V associated with S 23 shear loading, the matrix 
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I 



J 




IV 


V 


xz 


^xz 


IV 


V 


yz 


^yz 



2x2 



This may be written as 



,IV 

’xz 

0 



0 

1 



;V 

’yz 



(A-17) 



- 1 



(A-18) 



2x2 2x2 2x2 



From Eqs. (A-6) and, (A-18) the microstresses, {cr}, may be 

related to any specified shear loading as 

{a} =[t] {S} (A-19) 

2x1 2x2 2x2 2x1 

The derivation of Eqs. (A-16) and (A-19) is independent 
of the load state. For the particular case of plane macro- 
scopic stress treated in this investigation the macroscopic 
stress vectors become 

{S}*^ = <S^ 0 S2> 

(A-20) 

Equations (A-16) and (A-19) may be evaluated for the point 
by point microstresses associated with a plane macroscopic 
load condition as specified by Eq. (A-20) . 

Von Mises yield criterion is assumed to govern initial 
yield. This requires that 1 

+ T ['^= 






J = a ^ + o + 



xz 



yz ''xy 

- Vy - °y°z - = k 



1. + CT ^ + O ^ 

3 1- X y Z 



(A-21) 



Y is the elastic limit of the material for the tensile 
o 

stress state. 
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Equations (A-16) , (A-19 ) , (A-20 ) , and (A-21) may be 

combined to form a yield equation for plain macroscopic 
loading in terms of the macrostresses and the point depen- 
dent components of the A matrix. This equation is 



For given values of S and S„_ Eq. (A-22) may be solved 

X X z 

for the value of S necessary to cause yield. Each point 
within the fundamental block will give a different value of 
dependent on its microstress state. Each combination of 
and S^„ yields two values of S„ at each point, a maximum 
and a minimum. The minimum of the maximum values and the 
maximum of the minimum values determine two yield points for 
the given values of and Continuing in this manner 

the initial yield surface for a composite subjected to a 
macroscopic plane stress state is calculated. 





+ A|^) - I = 0 



(A-22) 
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APPENDIX B: CONSTITUTIVE EQUATIONS 

FOR A HOOKEAN FIBROUS COMPOSITE MATERIAL 



The well known equations of elasticity relating stress 
components and strain components in terms of two independent 
elastic constants represents a special case of generalized 
Hooke's Law for isotropic materials. Here we derive the 
Hopkean stress-strain relations for a lamina of an anisotropic 
fiber reinforced composite material. 

The generalized Hookean stress-strain law -relating stress 
and strain may be written in matrix form as 



e 




a , 


a ^ 


a 


a, , 


a. 


a. , 




S 


X 




11 


12 


13 


14 


15 


16 




X 






^21 


^22 


^23 


^24 


^25 


^26 




"y 






a 


a ^ 


a ^ 




a 






s 






31 


32 


33 


34 


35 


36 




z 


p 


= 


a , , 


a. , 


a . 










s 






41 


42 


43 


44 


4 5 


4 6 




xy 


p 




a 


a 


a 


a 


a 


a ^ 




S 


yz 




51 


52 


53 


54 


55 


56 




yz 


e 




a 


a 


a 


a 


a 


a 




S 


xz 




_ 61 


62 


63 


64 


65 


66_ 




xz 



In a more compact notation this may be written 



{e} = [aJ {s} 

6x1 6x6 6x1 



(B-2) 



The terms are engineering strain components and the S^j 

are stress components. 

Thus the generalized stress-strain relations contain 36 
unknown Hookean constants, 

a^j (i, j = 1,2,3,4,5,6) (B-3) 

Consideration of the strain-energy density function |27] 
establishes symmetry of the matrix of Hookean constants 



^ij ^ji 



(B-4) 
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and reduces the maximum number of independent elastic 
constants to 21 . Further reductions in the number of 
elastic constants depends on the number of planes of 
symmetry present in the material. 

When the elastic constants at a point have the same 
values for every pair of co-ordinate systems which are 
mirror images of each other dn a certain plane that plane 
is called a plane of elastic symmetry. This means that if 
a certain plane, say the (xj^,X2) plane in co-ordinate 
system (xj^,X2fX2) is a plane of elastic symmetry then the 
reflection of this system, (xj^,X2,X2) , where 




X3 = -x^ (B- 5 ) 

must have the same elastic constants. If you could actually 
"see" the elastic properties this means that the properties 
would appear identical at every point regardless of which 
side of the plane you were viewing. Alternatively Sines [2^ 
suggests that if you could actually see these properties 
they would look the same even if the observer stood on his 
head, that is rotated 180 ° about his line of sight. This 
invariance is utilized to evaluate the elastic constants. 

Consider the reference system (Xj^,X2/X2) of Fig. (B- 1 ) . 
If the plane (X2^,X2) is a plane of elastic symmetry a 
reflection through the x^ axis may be carried out without 
effecting the elastic constants. This new reference system 
is designated by the primed axes of Fig. (B- 1 ) . 
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Figure (B-1) . Reflection of a Plane of Syirunetry (x]^,X 2 ). 

The general equations governing the transformation of 
axes are 

^x " ^I'ex+^l'^^y+^l^ez+^l^'l^xy+^'iniEyz+ni^-ie^x 
^y = ^ 2 ' V2"^z+^2'^2^xy+"‘2^2£yz+^2^2ezx 
^z = V^3"^z+^3'^3^xy‘^“3^3^yz+'^3^3ezx 

^xy = 2\^2^ +2mim2^y +2nin2^z+ ( V2+"'l^2) ^xy 

+ (m3^n2+nim2) £y2+ (^1^2+^102) ^zx (B-6) 

C •* 

'•y2 = 2A2Jl3e^+2m2in3ey+2n2n2e2‘*' (^2^3''’^2^3^ ^xy 
+ (m2n3+n2m3)£y2+(^^2S + ^2^3)^zx 

^ZX = 2 I 13 I 11 ex+2m2inj^ey+2n3n2_e2+ (^'3n'i+^-ini3) e^y 

+ (m3n2^+n3mj^) £^2+ (^^S^l+^S^ii) Ezx 
for strain components and 
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Sy “ ^2^x'*‘^2^y‘^^2^z'^^^2^2^xy'^^^2^2^yz'*"^^2^2^zx 
^z ^ ■<^3Sx+n'3Sy+’^3S2'*'^^3"'3^xy+2m3n3Sy2+2n3£3S2x 
^xy = ^1^2Sx+"’l"*22y+'^l'^2^z+( V2+"'l^2)Sxy 

+ (min2+nini2)Sy2+(niJl2+5-1^2)2zx (B-7) 

Syz = ^2^3Sx+™2™3Sy+"'2^3Sz+(^2"'3+"'2^3)Sxy 
+ (m2n3+n2m3)Sy2+(n2J!'3+5<2n3)Szx 

^zx = A 35 ,iSx+m 3 inj^Sy+n 3 niS 2 +( 5-3111]^ +m 35, 

+ (ni3n;L+^3mi)Sy2+(^3Jil+^3ni)S2x 
for stress components. The terms 

1^, mj^, n^ (i= 1,2,3) (B-8) 

are direction cosines defined by the following table: 





1 — 1 


2 


3 


1' 




mi 


^1 


2 ' 


I 2 


m 2 


^2 


3 ' 


I 3 


^3 


^3 



This indicates that 1^^ cos (x|,X]^)etc. For the reflection 
of Fig. (B-1) the direction cosines are: 





1 


2 


3 


1' 


1 


0 


0 


2' 


0 


1 


0 


U) 


0 


0 


1 — 1 
1 
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Using these direction cosines to evaluate the transformation 
relations Eqs. (B-6) and (B-7) yields 



s • 

X 


= Sx 


1 

^x 


= 


^X 


s„' 


= s 








Y 


Y 


Y 




Y 


s ' 


= s 


e ' 


— 








z 
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s ' 


= -s 


e ’ 






xy_ 


xy 


xy 




xy 


S ' 


=-S 


e 


— 




yz 


yz 


yz 




yz 


s • 


= S 


1 


_ 




zx 


^zx 


ZX 




zx 



The stress-strain relation in the new co-ordinate system 
may be written 

1 I , I , T T 

<ejj Cy £2 ^xy ^yz ^zx ^ “ M *^^x ^y ^xy ”^yz ~^zx^ (B-10) 

6x6 

The elastic constants are then evaluated by equating the 
transformed strains to the original strains in accordance 
with Eq. (B-9) . Consider for example ej =£„♦ In detail 
this requires that 

aiiSx+ai2Sy+ai3S2-ai4Sxy-ai5Syz+ai6Sxz = 

^ll^x'*'^12^y'^^13^z'*'^l4^xy‘''^15^yz‘''^16^zx (B-11) 

The constants are invariant with transformation, requiring 



that 

“^15 ^yz ^ ^15 ^yz 
■^14 ^xy " ^14 ^xy 

For a non zero stress condition this can be true 

^15 = ®14 ' “ 

Carrying this out for the other relations of Eq. 
shows that 

^15"^^14=^24=^25=^34=^3 5=^46=^56^° 



(B-12) 
only if 

(B-13) 

(B-9) 

(B-14) 
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Thus a single plane of elastic symmetry reduces the 
number of independent elastic constants to 13. 

If a material has three orthogonal planes of symmetry 
a reflection of a second plane may be carried out without 
changing the Hookean coefficients. Reflecting the {X2r>^2^ 
plane through the axis results in the transformed axes 
of Fig. (B-2) . In this case the 



X 



1 




Figure (B-2). Reflection of a Plane of Symmetry (xo,x ) 

^ 3 

direction cosines are given by 



12 3 

1 ' 

2 ' 

3' 

The transformation relations Eqs. (B-6) and (B-7) give 



-1 


0 


0 


0 


1 


0 


0 


0 


1 
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(B-15) 



s' 
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e 
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x 


s' 
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S 
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z 
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-s 
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-e 


xy 




xy 


xy 




xy 
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— 
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yz 




yz 


yz 




yz 
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-s 


e ' 


= 




zx 




zx 


zx 




zx 



Using these relations it may be shown that 

^26 = ^16 = ^36 = ^45 = 0 ' 

The matrix of elastic constants may now be written 



Mo 

6x6 



^11 ^12 ^13 ° 

^22 ^23 ® 

^33 ° 

^44 

symmetric 



0 

0 

0 

0 



0 

0 

0 

0 

0 

^66 



(B-17) 



where the subscript o indicates that this matrix is defined 
for an orthotropic material. There are now nine independent 
elastic constants. 

Further reductions in the number of independent elastic 
constants depend upon the number of planes of isotropy, if 
any, present in the material. A plane of isotropy is one 
in which the elastic constants are invariant with respect 
to direction within the plane. 

Assume that a plane of isotropy does exist, the (Xj^,X 2 ) 
plane of Fig. (B-3) . This indicates that the Hookean 
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coefficients are independent of the orientation of the 
(x^,X 2 ) axes within this plane. An infinite number of 
orientations could be chosen within this plane to assist in 
evaluating the remaining unknovms . A 45 degree rotation, 
shown in Fig. (B-3) , was selected for ease of calculation. 




Figure (B-3). 45 Degree Rotation of Plane (xj^,X 2 ) 

The direction cosines for such a case are: 
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2 
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/2 


/2 
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/T 
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1 



Again using the transformation relations Eqs. (B-6) and 
(B-7) one obtains 



e ' = ^ 

X 


Ex + 
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^xy 


z^ = h 


^x + 


*5 


"y - 


*5 


^xy 



e' = e„ (B-18) 

X z 
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^xy ~ 



■ "/2 ^yz ^ 

Sx 



zx 



S' = - s + is + s 
X 2 X 7 y xy 

I + I Sy - S^y 



s' = s 

Z 2 

s' „ = i S -is 
xy 2 ^ 2 y 

S' =-^S +/is 

yz yz /2 zx 



,. _ 1 o . 1 



(B-18) 



^zx /2 ^yz ^zx 

Using these relations it may be shown that 



^55 


^66 


ail = 


^22 


^13 


^23 


It 


2(a 



(B-19) 



11 
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A material exhibiting one such plane of isotropy is 
called a transversely isotropic material. There are five 
remaining independent elastic constants and the matrix of 
elastic constants may be written 

0 
0 
0 



[* ..] ■ 



6x6 



^11 ^12 ^13 

^11 ^13 



'33 



2(aii aj^2^ 



0 

0 

0 

0 



symmetric 



55 



0 

0 

0 

0 

0 

‘55 



(B-20) 
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where the subscript TI signifies transversely isotropic. 

The stress-strain relations for this case may be written 

{e} = {S} (B-21) 

6x1 6x6 6x1 

One of the assumptions of this study is that a composite 
material may be considered to be homogeneous in a macro- 
scopic sense. This assumption allows the composite to be 
considered as a transversely isotropic material, and Eqs. 
(B-21) are the constitutive equations for a composite 
material. Introducing the engineering constants these 
equations become 

e = _L (Sj, - VS) - ^ S 

X ^ X i y Z 



"T 

e = — (S„ - v„s ) - 

y ' y T x' 

Vx 1 

^z = “— (S + S)+ — 

^ ^ X Y 



e = 1 c 

xy ^ ^xy 

£ = 1 C 

yz ^ yz 



(B-22) 



£ = 1 g 

zx ^ ^xz 



^T ■ 



2(1 + V^) 



E^ and are the Young's modulus and shear modulus in 
the transverse direction, in the plane of isotropy. E and 

Ij 

G are corresponding values in the longitudinal direction, 

Li 

perpendicular to the plane of isotropy. Poisson's ratio 
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v,p signifies the contraction in the plane of isotropy due 
to tension in the same plane. Poisson's ratio v represents 

1j 

the contraction in the plane of isotropy due to tension in 
the longitudinal direction. The transverse shear modulus 
is a dependent quantity. 
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APPENDIX C; YIELD LOCI 



The figures included in this section represent the 
yield loci that were generated during the course of this 
investigation. Each plot is labelled with the constituent 
materials, filament volume fraction, v^, and S and S 
scales. The outermost locus for all cases corresponds 
to Sjj 2 = 0* The increment between successive loci is 

also noted on each plot. 

Those plots that were used in the parametric analysis 
are labelled in terms of the parameter being investigated. 
All other properties are held constant. 

Note that in some cases the and S scales are dif- 

X z 

ferent from one another. This occurs when the matrix 
materials are quite soft in comparison to the filament 
material. The transverse strength is so much less than 
the longitudinal strength that the surface plots nearly 
as a straight line if the same scales are used for both 
axes . 
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(C-1) . Matrix Yield Strength Variation, Y„ = 3x10^ psi 
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(C-2) . Matrix Yield Strength Variation, Y = 2.3x10"* psi 
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(C-3) . Filament Stiffness Variation, = 30x10® psi 
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(C-4) . Filament Stiffness Variation, = 40x10® psi 
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(C-5) . Filament Stiffness Variation, = 75x10® psi 
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(C-6)., Boron/6061 Aluminum 
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(C-7) . Boron/6061 Aluminum 
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(C-8) . Boron/6061 Aluminxam 
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(C-9) . Boron/6061 Aluminum, Elliptical 
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Figure (C-10) . Boron/6061 Aluminum, Elliptical 
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(C-ll) . Boron/NARI4CO Epoxy 
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(C-12) . E-glass/Epoxy 
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(C-13) . E-glass/Epoxy 
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(C-14) . E-glass/Epoxy , 
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(C-15) . E-glass/Epoxy 
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(C-16) . E-glass/Epoxy 
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(C-17) . MODMOR I (graphite) /NARMCO Epoxy, 
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(C-18) . MODMOR I (graphite) /NARMCO Epoxy 
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(C-19) . MODMOR I (graphite) /NARMCO Epoxy 
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(C-20) . MODMOR I (graph! te) /NARMCO Epoxy, 



1 




(I 

> 



0 

u 

CP 

•H 



152 



(C-21) . MODMOR I (graphite) /NARMCO Epoxy, 
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(C-22) . MODMOR II (graphite) /NARMCO Epoxy, 
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(C-23) . MODMOR II (graphite) /NARMCO Epoxy 
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(C-24) . MODMOR II (graphite) /NARMCO Epoxy, 
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(C-25) . MODMOR II (graphite) /4617 Epoxy 
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(C-26) . MODMOR II (graphite) /4617 Epoxy, 
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(C-27) . MODMOR II (graphite) /4617 Epoxy, 
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(C-28). THORNEL 25 (graphite) /Epoxy 
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(C-29) . THORNEL 40 (graphite) /Epoxy 



APPENDIX D; PROGRA24 OVERLAY STRUCTURE 



The overlay root-segment structure is shown in Fig. 
(D-1) where each rectangle represents a subroutine. The 
subroutine functions are: 

MAIN -Remains in core and controls calling sequence 
during execution. 

BANDEC - Solves banded matrix problems. Used in MAINl 
and MAIN2. 

MULI -Multiplies a symmetric banded matrix with a 
vector. Used in MAINl and MAIN2. 

MAINl - Calculates microstresses and macrostresses for 
elastic loading of load cases I, II, and III. 

DAVES -Computes element stiffness and assembles system 
stiffness. 

CALTAU -Calculates element stresses and strains. 

REVGNK - Revises system stiffness to accommodate 
boundary conditions and applied loads. 

MAIN2 - Calculates microstresses and macrostresses for 
longitudinal shear elastic loading, load case IV. 

MAIN3 - Calculates Initial Yield Locus. 

MATINV - Matrix inversion with accompanying solution 
of linear equations. 
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Figure (D-1) . Overlay Root-Segment Structure 
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WRITE( 6,10 8) (I X,TCXX(1,IX) , TCYY ( 1 , 1 X ) , TCXY (1 , IX ) , TC ZZ ( 1 ♦ IX ) 
1IX=1.NUMSL ) 

WRITE (6,110) B(1,1),B{2,1),B(3,1) 
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OBTAIN THE B INVERSE MATRIX 

CALL MATINV( B» 3, S» 3 , 3 » IP I VOT , I MDS X tPI VOT*DET£RM) 
WRIT£(6,150) ((B(I,J)tJ=lt3),I=l,3) 
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